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ABSTRACT 


This fifth technical report, also the final report, covers studies 
on the effect of Impurities and defects on the performance of silicon 
solar cells which were not reported In the previous four technical 
rc:portS4 It describes a theoretical study of the effect of defects 
across the back-surface-field junction on the performance of hlgh- 
efflciency and thin solar cells, using a cfeueAopcd-p^A/imeXeA. device 
model for the three-dimensional defects. It shows that vevy signi- 
ficant degradation of open-circuit voltage can occur even if there 
are only a few defects distributed In the bulk of the solar cell. 

Two new features In the thickness dependences of the fill factor 
and efficiency in Impurity-doped back-surface-field solar cells 
are discovered in the exact numerical solution which are associated 
with the high injection level effect In the base and not predicted 
by the low-level analytical theory. What are believed to be the most 
accurate recombination parameters at the T1 center to date are also 
given and a new theory Is developed which Is capable of distinguishing 
an acceptor-llke deep level from a donor-llke deep level using the 
measured values of the thermal emission and capture cross sections. 
Using the measured emission-capture cross section ratio, this theory 
can also provide information concerning the magnitude of the lattice 
distortion around an Impurity atom before and after the capture or 
emission of an electron or a hole at the impurity center. 
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I . INTRODUCTION 


The objective of this progrea Is to determine the effects of Impurities 
and defects on the performance and permanence of silicon solar cells. It 
Includes theoretical (computer-aided-design techniques) and experimental 
studies of the effects of impurities on the electrical properties of silicon 
material Intentionally doped with specific impurity elements, and the effects 
of these impurities on the Impurity related energy level positions, the 
concentration of these energy levels and the recomblnatlon-generatlon-trapping 
rates of electrons and holes at these energy levels. 

This technical report is the fifth of five annual reports 
and it is also the final report of this contract. Studies performed 
prior to the fourth technical report have been reported In the first four 
technical reports and the results have been submitted for publication or 
published in engineering Journals. This report contains studies not reported 
in Che previous four reports. Reference to the work of the earlier technical 
reports are made from time to time but no attempt Is made to make a detailed 
review of these previous t^rork since they are documented In the four previous 
reports and in open literature. 

The first two chapters of this report (Chapter II and III) contain the 
results of a theoretical study of the effects of defects across the back- 
surface-field junction of hlgh-efflclency thin-base silicon solar cells doped 
with recombination Impurity elements. Since the geometry of the defect and 
Its effect on the cell are three-dimensional which defy exact numerical or 
exact analytical solution, a departure from the exact transmission line method 
employed In the previous four reports has been undertaken. A dtvctopeji 
peAAjn(iX.eA device model Is employed around t(ie defective regions of the cell and 
one-dlmenslonal low-level diffusion solutions of the dark current-voltage 


charactarlccics are applied to the regions of cell of this model. A two-regloo 
model Is developed for defects on the edge of the cell and a three-reglon model 
la developed for defects which are distributed In the bulk. Chapter II gives 
the results of the edge defects while chapter III ^ivea the results of the bulk 
defects. It is shown In Chapter III that volume distributed bulk defects 
across the back«>surface-fleld Junction can seriously degrade or reduce the 
open circuit voltage of a high-efflclency thin cell even only a few defects 
are present. The edge defects have less degrading effect If the area Is large 
or the area-to~clrcumfereace ratio Is large. However, performance degradation 
due to bulk defects Is relatively Insensitive to the area of the cell and 
depends only on the areal density of the defect. 

The third chapter of this report (Chapter iV) contains a short study 
of the effect of thickness on the solar cell performance parameters, the short- 
circuit current, the open-circuit voltage, the fill factor and the efficiency. 
The results were computed from the numerical solution of the Shockley Equations 
for one-dimensional cell structures wlUt and without back-surface field, using 
the tranmlsslon line circuit model technique. Hew features in the thickness 
dependence of the fill factor and the efficiency are presented for chick and 
very thin cells which are related to high Injection level and conductivity 
modulation conditions In the base region. These features were not observed In 
the thickness dependence curves given In Hovel's book which used the low-level 
analytical solution. 

Ttie fourth chapter of this report (Chapter V) contains a review of the 
literature data on the recorahlnatlon-generatlon-trapping parameters of electrons 
and lAoles at the T1 center In silicon and also gives the latest and we feel to 
be the most reliable data on the majority carrier thermal emission and capture 
rates and thermal activation energies at the two Tl levels. A new theoretical 
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development Is also presented to show how we can deduce the charge states 
or the acceptor-donor nature of the deep level using the measured nonequilibrium 
emission cross section and the measured equilibrium capture cross section of 
the majority carriers. It Is shown that this emission-capture cross section 
ratio can also provide information on the lattice relaxation around the impurity 
atom and provide the configuration entropy In a situation where the emission 
cross section Is nearly equal to Its equilibrium value. 

Plans are made to publish these four chapters as separate papers In 
open literature so that they are written as Independent chapters with their 
own literature references and equation and figure numbering systan. 


II. REDUCTION OF SOLAR CELL EFFICIENCY BY EDGE DEFECTS 


ACROSS THE BACK-SURFACE-FIELD JUNCTION 
A DEVELOPED PERIMETER MODEL ^ 

1. INTRODUCTION 

Back-surface-field solar cell contains a high-low junction near the 
back surface. The potential barrier of the space-charge layer of this high- 
low junction prevents the emit ter-j unction- Injected minority carrier from 
reaching the ohmic contact on the back surface. This greatly reduces the 
total base recombination current and gives a substantial Increase of the 
open-circuit voltage over a comparable cell which has no back-sur face-field 
high-low junction [1]. A typical result obtained using the computer-aided- 
design algorithm based on the transmlsslon-llne-clrcult model Is summarized 
In Table 1 [2]. It shows that the presence of the back-surface-fleld 


TABLE 1 EFFECT OF BACK- SURFACE-FIELD ON SILICON SOLAR CELL 
PERFORMANCE AT ONE AMI SUN AND 297*K. 


Computer 

Cell 

Impurity 

Profiles® 

VOC 

JSC 

FF 

EFF 

Run No* 

Type 

Emitter 

BSF 

mV 

mA/cm^ 


% 

RUN555 

P+/N 

exp(-x2) 

none 

554.11 

32.456 

0.8028 

16.236 

RUN335 

P+/N/N+ 

exp(-x2) 

exp(-X^) 

I 

672.32 

33.954 

0.7629 

19.587 


^Surface concentration®!. 25x10^°, bulk concentratlon=5xl0^Vm”^, thickness® 300pm. 


high-low junction, n/n+, on a p+/n junction silicon solar cell improves the 
open-circuit voltage by 672. 32-554. 11“118. 21 mV or a factor of 672.32/554.11 
=1.21. The AMI peak efficiency, EFF, is also improved by this factor, from 

16.236% to 19.587%. However, there is little increase of the short-circuit 

2 

current, only 1.498 mA/cm or 1.498/32/456*0.046 or 4.6%. This is expected 
since the cell thickness Is kept constant in this comparison so that the 
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# This chapter will be published in Solid-State Electronics, 1982 under the 
same title. 



total absorption or photocurrent density Is a constant. The fill factor Is 
decreased slightly when the back surface field Is added to the cell. This Is 
due to the high-level effect or higher minority carrier concentration In the 
BSF cell than the non-BSF cell because the low/high or n/n-<- junction reduces 
the minority carrier recombination rate, resulting In an Increase of the 
minority carrier concentration In the base region of the p+/n/n+ cell. 

In the above comparison, an Ideal Impurity diffusion profile, given by 
exp(-x /4Dt), Is assumed for boron In the diffusion p+ emitter and also for 

phosphorus In the dlffiised n4- back surface field layer. More realistic 

6 6 2/3 
diffusion profiles, such as exp(-x /Lj^) for phosphorus and [1 - ] 

for boron , which were used to study the effect of zinc recombination centers 

on silicon solar cell performance [2], do not change the conclusion that the 

BSF Increases the open-circuit voltage substantially. 

This dramatic improvement of the open-circuit voltage was first 
explained correctly in 1972 by Mandelkorn and Lamneck In series of experiments 
on silicon n+/p/p+ BSF cells whose top junction was made by phosphorus 
diffusion and whose BSF junction was made by alloying an evaporated layer of 
aluminum metal onto the back silicon surface [3]. A comprehensive theoretical 
analysis was then given by Godlewskl, Baraona and Brandhorst in 1973 [3]. 

They employed the one-dimensional model and the analytical solution is based 
on constant impurity doping in each layer, low level condition (minority 
carrier concentration much smaller than the majority carrier or dopant 
Impurity concentration) and finite surface recombination velocity at the 
metal-silicon interface on the back surface. In this n+/p/p+/m model for the 
experimental alloyed junction cells of Mandelkorn and Lamneck [3], the p+ 
layer contains two regions: the silicon regrowth layer which contains a high 
concentration of aluminum (7x10 to 1x10 Al/cm from Al-Si phase diagram) 


and an aluminum diffused silicon layer of about 1 urn thick from the 4-hour 
aluminum alloying at 800^0 . The m layer is the aluminum-silicon alloy which 
contains about 12 atomic percent of silicon. The surface recombination 
velocity is that at the p+/m Interface. It is evident that the Godlewskl- 
Baraona-Brandhorst model is only an approximation to the alloyed BSF junction 
solar cell. The GBB model is also an approximate model for an all-diffused 
BSF cell in which both the ni- emitter and p4* back surface field layers are 
obtained by Impurity diffusion. 

An earlier analysis of the back surface field was made by R. N. Hall 

in 1953 [4] who called a n/n+ junction in a p+/n/n+ diode a donor contact. 

The earlier applications were mainly for power rectifiers where the BSF 

layer would reduce the series contact resistance, increase the base layer 

conductivity modulation by getting into high injection condition at lower 

forward bias, and decrease reverse leakage current. 

In practice, the highest open-circuit voltage that has been obtained 

14 -3 

was from a silicon p+/n/n+ BSF cell made on 10 ohm-cm (5x10 cm ) float- 

zone silicon, reported by Fossum, Burgess and Nasby [5] in 1978. This cell 

2 

has an open-circuit voltage of 622 mV, short-circuit current of 34.3 mA/cm , 

2 

fill factor of 0.796 and efficiency of 16.8% under a 92.9 mW/cm simulated 
AMI solar illuminator. This experimental open-circuit voltage of 622 mV is 
substantially smaller than that predicted by the theory, 672 mV, shown in 
Table 1, Neither Interband Auger recombination nor energy gap narrowing in 
the heavily doped emitter seem to be able to account for the entire 50 mV 
difference between theory and experiments. 

For production n+/p/p+ back-surface-field silicon solar cells of high 
efficiency (14 to 16%), the open circuit voltage rarely exceeds 600 mV. For 
some of the earlier production cells, the back surface field p/p+ junction is 
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formed by aluminum alloying, similar to the method employed by Mandelkom 
and Lamneck [3]. More recent BSF silicon cells are produced by diffusion. 

A number of theoretical analyses has been advanced to account for the 
small experimental value of the open circuit voltage [6,7]. The general 
consensus Is that In these hlgh-efflclency silicon solar cells made on the 
hlgh-reslstlvlty base material (10 ohm-cm), base recombination rather than 
emitter recombination Is the limiting factor. The dominance of base recom- 
bination would place further emphasis on the Importance of employing the 
back-surf ace-fleld low-high junction to block the minority carriers from 
reaching the hlgh-recomblnatlon-rate ohmic contact on the back surface of 
silicon. This Importance of the BSF for the same reason Is also evident for 
thln-fllm cells In space applications due to their power/welght ratio 
advantage and for terrestrial applications due to chelr potential low 
manufacturing costs. 

The analyses made for the back-surface-field Junction In the past have 
focused on an areally uniform BSF high-low Junction. In the most general 
analysis made so far that has an explicit analytical solution, given by 
Godlewskl, Baraona and Brandhorst [1], the ohmic contact on the back surface 
is represented by a surface recombination velocity and the p+ layer of the 
p/p+ junction is assumed to have a finite minority carrier lifetime. These 
realistic assumptions made the lot.-hlgh Junction not a perfect minor Ity- 
carrler-ref lectlng potential barrier. However, the potential barrier of the 
low-high Junction is usually several tenths of electron-volt or many kT/q. 
Such a 'high' potential barrier would reduce the effect of both the bulk 
recombination In the p+ layer as well as the surface recombination on the 
back surface. Thus, a well-designed and fabricated low-high Junction would 
usually shield the minority carriers almost completely from these regions 
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of high recombination rates, the p+ layer and the p4>/m Interface. 

This ideal areally homogeneous model for the BSF junction overlooks 
some very Important defective conditions which Invariably occur In production 
or even In experimental solar cells, due to cell geometry, material defects, 
as well as inhomogeneities. These defects are divided Into two groups according 
to their locations and are Illustrated In Fig. 1: (A) perimeter or edge defects 
and (B) bulk defects. The origin and causes of these defects are given below. 

(A) Perimeter or Edge Defects 

(1) Overflow of solder 

(2) Imperfect perimeter surface etch 

(3) Poor encapsulation 

(B) Bulk Defects 

(1) Poor alloying of the BSF junction 

(2) Poor diffused BSF Junction due to diffusion pipes 

(3) Random back contact metal penetration 

(4) Impurity segregation and cluster 

These defects are electrical short circuits or resistive shunts across 
BSF junction. They can greatly reduce the effectiveness of the BSF potential 
barrier to shield the minority carriers from the recombination centers in the 
p+ region and on the p+/m interface. 

These defects are even more detrimental to the cell performance if 
they appear across the front junction. However, their degradation influence on 
the front junction is well known and manufacturing procedures are developed to 
minimize their appearance across the front junction. Their influence on the 
performance of BSF junction have not been generally recognized nor analyzed. 

In this paper, we shall analyze the effects of these electric defects 
across the BSF junction. We shall show that they can significantly reduce the 
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open-circuit voltaget hence the efficiency. 

The device models used In this analysis are analytical models rather 
than the exact computer models since the localized defect problem Is Inherently 
a three-dimensional problem which cannot be solved exactly and numerically at 
present due to computer size and speed limitations. Although the models to be 
used are approximate one-dlmenslonal analytical models, the results should 
provide accurate qualitative understanding and good quantitative estimate of 
the effects of defects on the solar-cell open-circuit voltage and performance. 

Of the two locations of the defects, edge and bulk, the edge defects 
will be analyzed In this paper. Both a perfect short-circuit par.ch around the 
entire perimeter and an Imperfect patch characterized by a finite surface 
recombination velocity will be analyzed. The second group of defects, present 
In the bulk, will be analyzed In a subsequent paper [8]. 
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2. ENTIRE PERIMETER SHORT CIRCUITED - DEVICE MODEL AND ANALYSES 

The back surface field la nomally a very thin layer of about one 
micrometer thick fabricated either by aluminum alloying In earlier n4>/p/p^ 
silicon cells or by diffusion In recent p^/n/n4- and n4'/p/p+ silicon cells. 

The rather thin back surface-field layer Is mainly from two reasons. It Is 
fabricated at a low tonprature (alloying at 800*C or diffusion at 900 to 
1000"C) to preserve the minority carrier lifetime In the base region. This 
low temperature process gives a thin layer. The second reason Is that the 
effectiveness of the low-high Junction as a potential barrier for the 
emitter-injected minority carrier is the best tdien the low-high Junction 
transition occurs In a very small distance or Is very abrupt. This would 
mean that the n/n4- low-high Junction should have a high Impurity concen- 
tration In the n4- layer as well as a very thin space-charge layer to give a 
large and abrupt potential barrier across the n/r& Junction. 

In practice, the perimeter of this thin back-surface-field layer can 
be short-circuited by the alloy used to solder the cell to a heat sink. Such 
a perimeter short-circuit can significantly reduce the open-circuit voltage 
since these shorts will substantially Increase the minority carrier recombina- 
tion current In the perimeter region. In this section, we shall treat the 
perimeter patch by a short-circuit so the perimeter of the cell is degraded 
from a p+/n/n+ diode to a p+/n diode. In a later section, a finite surface- 
recombination velocity will be used to characterize the perimeter shunt. 

This is a two-dimensional problem if the cell Is perfectly circular. 

It is a three-dimensional problem If the cell has a square, rectangular or 
other geometries. In the low-level and content bede-dopant-impurlty 
concentration case, the analytical solutions In series forms can be obtained 
for the mlnolrty-carrler diffusion equation In the base region. However, the 
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analytical solutions cannot be extended to nor used for a cell with a random 
distribution of patchy defects across the perimeter or the bulk of its back- 
surface-field low-high Junction. Thus* we shall use an approximate model by 
dividing a cell of arbitrary shape into two regions, the perimeter region and 
the central region and use the one-dimensional analytical solutions for these 
two regions. 

Two examples of this two-region model are Illustrated in Figure 2. 

One is the square cell and the other is the circular cell. The mathonatical 
models of the two regions are as follows. The central region (Region I) is 
assumed to have the original p+/n/n+ BSF cell characteristics. The perimeter 
region (Region II) Is assumed to experience the full Influence of the short- 
circuit patch around the entire perimeter of the low-high Junction due to 
soldering shorts or mechanical damage. This Is the extreme or worst case of a 
defective perimeter cell. The Intermediate case is modeled by a finite 
surface recombination velocity across the perimeter, analyzed in a later 
section. 

From Fig. 2(b), It is evident that the perimeter region has a triangular 
cross section. The width of this triangular belt Is assumed to be given by 
as Indicated in Fig. 2(b). It is the width of the p/n Junction of 
this perimeter region. Here, Is a few times of the base minority carrier 
diffusion length, where is the minority carrier diffusion 

length whose dlffusivity Is and lifetime n is a number of the order of 
one and In the numerical solutions to be presented graphically, n Is varied 
from 1 to 4 to show that Its value does not have a large effect on the 
degradation of the open-circuit voltage. would become the amblpolar 
diffusion length If we are In high Injection level. In Fig. 2(b) Is the 
thickness of the base region and It is approximately the cell thlukuess since 
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the diffused emitter and baclc-surface-fleld layers are usually very thin 
compared with the cell thickness. 

The significance of L^, which Is the length of the side of the 
triangular cross section of the perimeter region as shown In Fig. 2(b), Is 
as follows. All the minority carriers injected by the p/n Junction will be 
Influenced by the short-circuit at the defective perimeter If they are injected 
within a distance of from the short-circuit. We call the distance of 
influence or the range of the defect. A one-dlf fusion-length model is one 
with n«l or Lj 

perimeter short-circuit since when the cell is thick or X_ > L_“L_, we will 
have no edge short-‘circuit effec , even If the low-high Junction around the cell 
perimeter is completely short-circuited. Thus* the one-dlf fusion-length model 
will provide a lower estimate of the degradation effect on the open-circuit 
voltage due to perimeter short-circuits. 

The current-density equations for the minority carriers injected by 
the emitter Junction in regions I and 11* which are now treated as one-dimensional 
diodes* are the well-known Shockley diode equation at low injection level* 
modified by the effect of the ohmic (infinite surface recombination) short- 
circuit in region II or the back surface-field in region I. These are given by 

Jf - q(DgpQ/Lg)tanh(Xg/Ljj) (1) 

and 

Jjl * q(DgpQ/Lp)ctnh(X/Lg). (2) 

The synods here follow the conventional usuage. PQ"Pg(exp(qV/kT)-l] . 
q Is the magnitude of the electron charge, k Is the Boltzmann constant. 

T is the cell temperature, P_ is the minority carrier concentration 
in the uniformly doped base* and V is the terminal voltage or the 
difference between the quasi-Fermi or electrochemical potentials of 
boles and electrons across the emitter-junction space-charge Layer. 


«nLg*Lg. Using n«l would underestimate the effect of the 
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JL^ Is tbs bsss Isyer thlckosss. X is ths dlstsncs bstvssn sn slsnsntsl junction 
area in Ragion II and the nearest short-circuit point at the perineter as shoim 
on the right cross-sectional picture in Fig. 2(b). Thus* Region II has an 
effective variable base layer thickness, X(y) , t^ere y is a direction in the 
plane of the p/n junction. 

current density of Region I, given by Equation (1), assumes that the 
low-high junction is perfectly minority-carrier-reflecting, that is, it has an 
infinite potential barrier height so that no minority carriers can cross this 
low-high junction boundary. This assumption is an excellent approximation for 
all practical cases of well designed and carefully fabricated cells. Neither 
Fermi-Dirac statistics nor energy-gap narrowing in the diffused or alloyed pv 
back-surface-field layer of a nf-/p/p4> diode, or in the nf diffused BSF layer 
of p^/n/nf diode, can change the numerical results given by Eq. (1) significantly. 
This can be demonstrated numerically by evaluating the more complete equation 
of given by Godlewski, Baraona and Brandhorst [1] and assuming (1) a perfect 
ohmic contact (Infinite surface recombination) on the back surface of the cell 
and (11) as much as 100 mV reduction of Che energy gap in the n^ or pf region 
of the BSF layer. 

The total diode current from the two regions are obtained by integrating 
the current densities. Equations (1) and (2), over the cross-sectional areas of 
the two regions. With the aid of Figs. 2(a) and 2(b), they are given by 

ll - q(DgPQ/Ljj)tanh(Xg/L^)A. (3) 

and 

^2 * j ctnb(X/L^)dx (4) 

0 

The symbols used here were defined after Equations (1) and (2). 

The current in Region I, given by Equation (3), is just the current 
density, J^, times the area of Region I, A^. 
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Th« current In Region II is given by Equs.tlon (4) where is the 
effective length of the circtfinfcrence of the perimeter region which is defined 
by the ares of Region II after Equation (4) is avsluated. For simple geometries 
such as the square and the circular cells* it can be obtained geometrically 
and is 

Zij - HZ - >q^) <5) 

for a square cell and 

Zjj - (D - - Xj) (6) 

for a circular cell. Here, Z is the edge dimension of the square cell and D 
ia the diameter of the circular cell. 

The integral in Equation (4) can be evaliuited numerically. Since the 
model is approximate and an oae-dimenslonal approximation was already employed 
in the current densities given by Equations (1) and (2), the integral in 
Equation (4) is also evaluated approximately to give an analytical solution. 

To this aim, the triangular cross section of Region II is un^otdcd or dtveZopCd 
into a trapezoidal cross section shown in Fig. 2(c). This is reminiscent of the 
developed helix model used in the early theory of helix traveling wave tubes (9J. 
The trapezoidal or linear transformation for region II is given by 

X(y) « X - + I(Lj-Xg)/»^Lj-X2]y ( 7 ) 

which is used to evaluate the integral in Equation (4) in order to give an 

analytical solution. In the appendix, it is shown that the linear approximation 

gives numerical solutions which differ frexn the exact evaluation of the Integral 

2 2 2 

using the quadratic transforaatlon, X « X^ -f y , by less than 14Z. 

After regrouping terms and taking out the voltage dependent factor, 
(exp(qV/kT) - IJ, the total current flowing in the these two regions are given by 


and 


I 2 - J2A2lexp(qV/kT) - 1] 


(9) 


Here, the current density coefficients are defined by 


•'i ■ 


Jj - tqDjPj/(lj-X,)J log^lsll>httj/lj)/aIilh(Xj/Lj)l, 


The current flowing In an Ideal cell without any edge short Is given by 

Iq - J^AQ[exp(qV/kT) - 1]. (12) 

The areas, A^, A^ and A^, are the total area and the area of regions I 
and II respectively. A^^A^+A^. For a square cell, these are given by 

Aq « (13) 


Aj^ - (Z - 


[4(Z - - XJ)J 


For a circular cell with diameter D, these areas are given by 
Aq = irD^/4 


* Tr(D - 


^[7T(D - ^ 


Note that the area ratio, circuit and square cells are identical* 

This is also true for the area ratio A^/Aq. 

The degradation of the open-*circuit voltage, -AV^^, can be ccmiputed by 
ccmparing the total current of a perfect cell, 1^, with the total current of a 
defective cell, The result is independent of the cell operating 

condition or the magnitude of the open-circuit voltage if the open-circuit 
voltage is large compared with kT/q or about 25 mV so that the approximation, 
[exp(qV/kT) - 1] - exp(qV/kT), is accurate. For a photocurrent or short-circuit 
current of from the solar Illumination, the defective cell gives 

II = Ii + I2 = 


and the perfect cell gives 




J^A- {exp[q(V-_/kT)] - 1} 


Dropplsg 1 compared with exp(qV^^/kT) and eliminating between Equations (19) 
and (20) » the reduction of the open-circuit voltage Is then given by 

- AVqp - (kT/q) log^KAj^/Ag) + (J 2 A 2 /JjAq)J (21) 

where and J 2 are given by Equations (10) and (11) and the areas are given by 
Equations (13) to (15) for the square cell and (16) to (18) for the circular or 
round cell. 

As stated after Equation (18) > the area ratios A^/A^ and si's 

Independent of whether the cell Is square or round. Thus, the degradation of 
the open-circuit voltage, given by Equation (21), Is the same for these two 
cell geometries using the developed perimeter model. However, we would expect 
degradation to have some dependences on the cell geometry In general but 
such differences are probably not very large for the perimeter or edge defects. 

3. EFFECT OF THICKNESS IN PERIMETER DEFECTIVE CELLS 

The analytical solution for the developed perimeter model just obtained 
provides a quantitative estimate of the Importance of a short-circuiting defect 
across the entire perimeter of the BSF junction on the open-circuit voltage 
and cell performance. When the cell Is thin compared with the diffusion length, 
Xj « Lg, the BSF layer gives a high compared with a cell without the BSF. 
Thus, we would expect the perimeter short-circuit to have a larger reduction of 
as the cell becomes thinner. From a simple geometrical consideration of 
Figure 2, It Is also evident that V „ reduction by edge defect Is larger If the 
cell area is smaller or if the edge dimension Z or the diameter D is not much 
larger than the base minority carrier diffusion length, i.e*, Z if or 
if These are physically obvious conditions under which the perimeter 

short-circuit will give a large reduction of the open-circuit voltage* They 
may also be illustrated numerically by the following examples* 



2 

Consider an experimental cell of 1 cm (Z-lcm or D-lcm) and L^>1000 ym 

^hlch is the order of that reported for the hlghest-V^g dlffused-slllcon cell 

[5]. Let us also assume a cell thickness of 500 ym, then we have Z/L «10, 

B 

Xg/Lj-0.5, Aj^/Aq-0.6836, A 2 /AQ-O. 3 I 64 and J2/Jj^-3.52 using (10) to (18). From 
(21), the computed open-circuit voltage reduction at 300”K is 

- - (kT/q)log^(1.797) - 0.58(kT/q) - 15.2 mV. 

If the cell thickness is reduced to about 50 ym, which is approximately 
the optimum thickness for a thin BSF cell to get the maximum AMI efficiency 
of about 17% [10], then Xg/Lg-0.05, Z/Lg=10, A^/Aq=0.6404, A2/Aq= 0.3596 and 
J2 /Jj^“ 63.19 using (10) to (18). From (21), the open-circuit voltage reduction 
would be increased to 

- * (kX/q)log^(22.72) » 3.15(kX/q) = 81.46 mV 

at 300"K. Xhis is quite a significant reduction or degradation of the open- 

circuit voltage by the presence of a defect. 

Xo provide a rapid estimate of the effect of perimeter short-circuit, the 

amount of open-circuit voltage degradation is computed from Equation (21) and 

presented In Figure 3 (solid curves with L /L =1.0). Xhe - AV are given In 

X B UC« 

10 -3 

mV at 297.15®K (n^=10 cm for Si) as a function of the normalized cell 

thickness, (from 0.01 to 1, 2 or 4), and using the square root of the cell 

area or the cell linear dimension as the constant parameter ox* D/L^ = 10 to 

100). It is evident that the - becomes more important in thinner cells. 

The different choices of the distance-of- influence are also investigated. 

The results of - L^/Lg=l (solid curves), 2 and 4 (broken curves) are 

compared in Figure 3. It is evident from the three upper curves for Z/Lg=10 and 

the two lower curves for Z/L =100 that the choice of L- have little effect in 

B X 

thin cells when Xg/L^«l. For thick cells or near Xg/L^-1, using L^/L^=2 would 
give better estimate since Lj./L^=l gives zero However, in thick cells, 

small anyway as expected physically. 
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AV^^(mV) 


THICKNESS Xb/L, 


The reduction of the open-circuit voltage as a function of 
the cell thickness of back-surface-field solar cells with 
edge defects across the BSF junction. The constant 
parameter is the cell diameter or cell edge ncrraalized to 
base minority carrier diffusion length. 
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IMPERFECT EDGE SHORT CIRCUITS 


The perfect short circuit across the entire perimeter of the low-high BSF 
junction, analyzed In the preceding section, can be extended to a more general 
situation In which the edge defect Is not a zero-resistance or perfect short 
circuit but has a finite surface recombination velocity, to be denoted by S^. 

The analyses Just made for the perfect short circuit can be easily extended to 
this case using the two-region, one-dlmenslonal analyses applied to the developed 
cell model. The solutions for the current density, J^, and total current, I^, 
are still given by Equations (1) or (10) and (3) or (8) respectively. 

The solution for the perimeter region. Region II, given previously by 
Equations (2), (4), (9) and (11), needs to be modified. With a finite surface 
recombination velocity, S^, the current density In the triangular region II is 
given by 

[sinh(X/L_) + (S_L_/D )cosh(X/L_)] 

J__ « q(Vn/^> ^ — <22) 

® ° ® tcosh(X/Ljj) + (SjjLg/Dg)sinh(X/Lg)] 

which replaces Equation (2) and reduces to Equation (2) if for a perfect 

short circuit across the perimeter of the BSF low-high junction. 

The ratio of the hyperbolic functions given in (22) replaces ctnh(X/Lg) 
of Equation (4) and It can be Integrated analytically since (22) is an exact 
differential if the linef.r transformation given by Equation (7) is used. The 
result can be put into the form of Equation (9) but the current density 
coefficient, J 2 , Is now given by 


J 


2 




cosh(Lj./Lg) + (SpLg/Dg)slhn(Lj./Lg) 

cosh(Xg/Lg) + (SjjLg/Dg)sinh(Xg/Lg) 


(23) 


The asymptotic solution of J 2 given above for Sp=<» reduces to Equation (11), the 
correct asymptotic behavior. However, Equation (23) gives incorrect J 2 when Sp=0 
since in this limit, there is no open-circuit reduction. In order to 
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give the correct asymptotic solution at S^"0, the cosh terms In Equation (23) 
are repalced by exp[ (Lj/Lg)tanh(X^/Lg)] and exp[(Xg/Lg)tanh(Xg/Lg)J so that 




exp[(Lj/Lg)tanh(Xg/Lg)] + (SjjLg/Dg)slnh(Lj/Lg) 

exp[(Xg/Lg)tanh(Xg/Lg)] + (SpI^/Dg)slnh(Xg/Lg) 


( 24 ) 


The reduction of the open-circuit voltage can now be computed from 
Equation (21) using and A 2 before and the new J 2 given by Equation (24). 

The normalized surface recombination velocity, S L /D , lies in the range of 

D D o 

2 3 2 4 2 

10 to 10 since in practice, D =1 to 10 cm /sec and L„=10 to 10 cm, ^lle 

ii D 

2 6 

Sjj==10 to 10 cm/sec. Thus, the diffusion velocity defined by lies 

3 4 

generally in the range of 10 to 10 cm/sec. Combining these with the extremes 
of S^, we have St./(D^/L^)= 10^/10^-10 ^ at the lower end and 10^/10^=10^ at the 
upper end. 

In Figure 4, the reduction of the open-circuit voltage is plotted as a 

function of cell thickness for a cell size of Z/L *10 and for surface 

B 

2 -1 

recombination velocities of S_L /D » «, 10 , 10, 1, 10 and 0. This family of 

D B B 

curves show that if the defect across the perimeter of the BSF junction is not 

a perfect short circuit, the open-circuit voltage reduction is decreased. As a 

2-1 3 

numerical example, let D =10 cm /sec, L =10 cm and S =10 cm/sec. Then, the 

B B 2 

2-13 

diffusion velocity is D /L =10 /lO =10 cm/sec so that S L /D =1.0. For a cell 

B B D B B 

of 100 ym thick, the open-circuit voltage is reduced by 32 mV by a imperfect 
perimeter defect compared with 60 mV if the defect is a perfect short circuit or 


2 2 

The above example has a rather small area (A^=Z = 1 cm ) and Z/Lg=10. 

For a larger cell, the open-circuit voltage degradation by defective perimeter 

will be smaller. However, the degradation can be decreased by a factor of two 

if the surface recombination velocity is decreased from a nearly perfect short- 

6 3 

circuit value of 10 cm/s to a medium value of 10 cm/s as indicated by the 
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THICKNESS Xb/Lb 


Fig. 4 The reduction of the open-circuit voltage as a function of 
the cell thickness of back-surface-field solar cells with 
imperfect edge defects across the BSF junction. The cell 
diameter or size is ten times the ba.T^: minority carrier 
diffusion length. The constant parameter is the effective 
recombination velocity of the edge defect, Sj^, normalized 
to the diffusion velocity of the minority carriers in the 
base, Dg/Lg. 





3 2 

two curves In Figure 4 with S^L^/D^ ■■ 1 and 10 (between 10 and «). 

In Figure 5, the reduction of the open-circuit voltage Is also 
presented as a function of the normalized surface recombination velocity 
with the cell thlcknes, as the constant parametar. It Is evident from 

this family of curves that open-circuit voltage degradation begins to saturate 
when the surface recombination velocity of the perimeter defect Is greater 
than ten times diffusion velocity, l.e., S_L_/D. > 10. Both this figure and 

Figure 4 show that the - Is less than (kl/q) or 26 mV when the cel] 

thickness is greater than about 0.35L- In a cell whose BSF Junction perimeter 
Is completely short-circuited. 

5. SUMMARY AND CONCLUSION 

The effects of electrical short circuits across the back-surface-field 
junction at the perimeter or the edge surface of a solar cell Is analyzed by 
applying the low-level one-dimensional analytical theory to a 'developed-cell' 
geometrical device model using the concept of distance-of-influence for the 
perimeter region of the cell. This perimeter or edge defect can cause a 
significant reduction of the open-circuit voltage and efficiency of a BSF cell. 
For small area, thin and large base-diffusion-length cells, the reduction of the 
open-circuit voltage and efficxency can be very large, so large that the edge 
defects can almost completely nullify the performance improvement offered by the 
back-surface-field junction. The results also Illustrate the Importance of 
edge or perimeter defects on the variation of the open-circuit voltage and 
efficiency In production lots of high-efficiency BSF solar cells. 
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Fig. 5 The reduction of the open-circuit voltage as a function of 

the defect surface recombination velocity of solar cells with 
imperfect edge defects across the BSF Junction. The cell 
diameter or size is ten times the base minority carrier 
diffusion length. The normalized cell thickness, X^/L^, is 
the constant parameter. 
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APPENDIX 


In this appendix* the accuracy of the developed perimeter approximation, 

used to evaluate the Integral given by Equation (4), Is checked against an 

exact numerical evaluation of this Integral. The exact Integral uses the 

circular transformation 
2 2 2 

- Xj + y^ (A.l) 

Instead of the linear transformation given by Equation (7). Thus, the exact 
and the approximate expressions for the Integral given by Equation (4) are 


and 


I. (appro) ■ i/(L+W)/(L-W) log [sinh(L)/sinh(W)] 
L e 

I-(exact) ■ I ctnh(»'9^4W^) dy 

Jq 


(A.2) 


(A. 3) 


where L»Lj/Lj and W»Xg/Lg. Equation (A.2) Is similar to where is 

the effective length of Region II given by the ( ] term in Equation (15) or 
(18) for the area A 2 . 

The percentage error, defined by 100[1 - l 2 (Appro)/l 2 (exact)], Is 
computed numerically as a functon of cell thickness with the dlstance-of- 
Influence, L_/L_, varied from 1 to 4. The largest error occurs in the thin 

1 D 

cell range where X^/Lg^O.Ol to 0.2. The error varies from 10% to 14%. As 

X-/L- Increases beyond 0.2 towards 1.0 or L /L , the error decreases rapidly 
D B IB 

towards zero. The results are graphed In Figure 6. 
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THICKNESS Xb/Lb 


Fig. 6 The percentage error of the open-circuit-voltage computed from 
an analytical approximation to the integral* plotted ao a 
function of the ce^’’ thickness. 
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III. REDUCTION OF SOLAR CELL EFFICIENCY BY BULK DEFECTS 
ACROSS THE BACK-SURFACE-FIELD JUNCTION^ 

1. mit(»)UCTION 

Back-surface-field aolar cells contain a high- low (r4*/o or p^/p) 
jiuictlon in the semiconductor near the back surface on which a metal- 
semiconductor ohmic contact is made. The potential barrier of the space- 
charge layer of this high^low Junction prevents the emitter-junction injected 
and photogenerated minority carriers from reaching the ohmic contact on the 
back surface. Without the back-surface field, the minority carriers would 
reach the ohmic contact on the back surface and recombine with majority carriers 
at very high recombination races. This would greatly increase the recombination 
current in a solar cell and decrease the efficiency and the open-circuit voltage. 
Without a back-surfsce field, the cell performance would decrease as the cell 
thickness is decreased since the ohmic high-recomb inst ion-rate back metal- 
semiconductor Interface would be closer to the front injecting and collecting 
p/n Junction. Thus, to achieve high efficiencies in solar cells, modern cell 
designs always contain a back-surface-field high- low Junction. 

The presence of a back-surface-field Junction can increase the cell 
performance substantially over a cell without the back-surface-field Junction. 

We shall call the cells without the back-surface-field Junction the bocfe- 
■iuA^acc^ohnic cells and use the acronym BSF for back-surface-field and BSO for 
back-surface-ohmlc cells in this paper. An illustration of performance 
Improvement by the presence of a BSF is given in Fig. 3.1 where efficiencies 
of the BSF cells type) and the BSO (n4-/p type) are given as a function 

of the base layer thickness, X^. These curves were computed by numerically 
solving the one-dimensional Shockley Equations using the transmission line 
circuit model technique (Ij for a low-level base minority-carrier diffusion 
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# This chapter will be published in the April 19S2 issue of the Journal of 
Applied Physics under the same title. 


length of 577 pm under one AMI (terrestrial) solar lllunlnatlon. Such a long 
diffusion length was tised to simulate state^f-the-art high-efficiency silicon 
solar cells [2]. 

^ig« 3.1 sho%ra that the BSF cells have significantly higher efficiencies 
than the BSO ce3.1s as the thickness decreases. At 100 urn cell thickness, the 
BSF cell efficiency is 16.3% - 12.3% - 4% higher than the BSO cell. 

The principal cell parameter that causes the lower efficiency In the 
BSO cells Is the r duction of the open-circuit voltage In the BSO cells due to 
high recombination rate at the back-surface ohmic contact. This Is illustrated 
In Fig. 3.2 for the two types of cells whose efficiencies were given In Fig. 3.1. 
In Fig. 3.2, the open-circuit voltage tinder one AMI sun is given as a function 
of the cell thickness for the BSF and BSO cells. As an example, the open-circuit 
voltage In the back-surface-field cell of 100 pm thickness Is 607 - 495 » 112 mV 
higher than the back-surf ace-ohmlc cell. 

The very large reduction of the cell performance due to the lack of a 
BSF high-low junction on the back surface just Illustrated suggests that random 
defects across the BSF junction can significantly reduce the performance of 
hlgh-efflclency cells and give rise to a large spread of cell characteristics 
In production lot quantities. The random defects that give rise to resistive 
shunts across the BSF low-high junction can be put Into two groups: those that 
appear across the perimeter of the low-high BSF junction or on the edge of the 
cell and those that appear In the bulk of the low-high BSF junction. These are 
graphically shown In Fig. 3.3. The principal causes of these edge and bulk 
defects are mainly from the many manufacturing steps during the fabrication of 
the cells and they are tabulated In the following table. Their presence 
will reduce the cell performance from that of the Ideal or perfect BSF cells, 
given by the upper curves In Fig. 3.1 and 3.2, to the shaded region below. 


- 32 - 


EFFICIENCY 




THICKNESS(iuM) 


AMI efficiencies of BSF (back-surface-field) n+/p/p+ and 
BSO (back-sur face-ohmic) n+/p silicon solar cells as a 
function of cell thickness with 577 ym base minority carrier 
diffusion length, computed by exact numerical solution of 
the one-dimensional Shockley Equations. Defective cells 
with defect across the BSF junction have efficiencies between 
the BSF and BSO cells. 







THICKNESS(iuM) 


Fig. 2 AMI open-circuit voltage of BSF (back-surface-field) n+/p/p+ 
and BSO (back-surface-ohmic) n-f/p silicon solar cells as a 
function of cell thickness with 577 ym base minority carrier 
diffusion length, computed by exact numerical solution of the 
one-dimensional Shockley Equations. Defective cells with 
defect across the BSF junction have open-circuit voltages 
between the BSF and BSO cells. 







(A) Edge Defects 

(1) Overflow of solder 

(2) Imperfect perimeter surface etch 

(3) Poor encapsulation 

(B) Bulk Defects 

(1) Poor alloying of the BSF junction 

(2) Poor diffused BSF junction due to diffusion pipes 

(3) Random back contact metal penetration 

(4) Impurity segregation and cluster 

These defects are even more detrimental to the cell performance if they 
appear across the front p/n junction. &>wever, their degradation influence on 
the front junction Is well known and manufacturing procedures are developed to 
minimize their presence across the front junction. Their Influence on the 
performance of BSF junction have not been generally recognized nor analyzed. 

In another paper, we have analyzed the effects on solar cell performance 
from edge defects across the BSF junction [3]. It is shown that these edge 
defects can significantly reduce the open-circuit voltage in thin hlgh-efficlency 
cells. The reduction decreases when the area of the cell Increases since the 
edge defects In a large area cell Is less effective as a hlgh-recond)lnation 
sink for the Injected minority carriers at the central portion of the cell 
which are f ar . away . from the edge of the cel'* In large area cells. However, 
the bulk defects across the BSF junction are distributed over the entire area 
of the cell and hence their performance-reduction Influences cannot be reduced 
by Increasing the area of the cell. as in the case of edge defects. Thus, one 
would expect a much higher performance reduction caused by bulk defects than 
by the edge defects across the BSF junctions of high-efflclency solar cells. 
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The degradation of the solar-cell performance due to a distribution of 
bulk defects across the back-surface-field Junction is analyzed in this paper. 
The device model employed is described in section 2 and the nathesiatical 
formulae are also given there. A detailed derivation of these formulae is 
given in the appendix. Families of curves are computed from these formulae 
and their physical significances are discussed in detail based on the three 
parameters which are used to characterize the defects: defect area, defect 
density and defect surface recombination velocity. The reduction of the open- 
circuit voltage due to the presence of the defect is presented as a function of 
defect area, density, cell thickness and defect surface recombination 
velocity. In section 3, the effects of short-circuiting defects are 
analyzed while in section 4, the effects of finite defect surface recombination 
velocity are described quantitatively. Numerical examples are given in each 
section to Illustrate the importance of the particular defect parameter using 
nisnbers that are similar to those of the state-of-the-art high-performance 
silicon solar cells. A summary of this theoretical study is given in 
section 5. 
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2 . 


DEVICE MODEL AND ANALYSES 


The edge and bulk defects across the BSF junction of a solar cell are 
three-dimensional defects. They have random sizes and are randomly distributed. 
Thus, the nonlinear Shockley Equations which describe the diffusion, drift, 
recombination and generation of electrons and holes in a semiconductor cannot 
be solved analytically or numerically for solar cells which contain these random 
defects. Exact numerical solutions could be obtained for two limiting cases: 

(1) there is a negligible number of defect and (2) there is a very high density 
of defects which completely short circuit the BSF Junction. These two ideal 
cases were computed and used as Initial Illustrations given in Figs. 3.1 and 3.2 
to show the importance of the defects on the cell performance. 

In order to obtain quantitative results for intermediate defect densities 
so that their effect on solar cell performfmce can be numerically estimated, 
approximate device models must be employed. A developed pM/iine^£A device model 
was constructed to analyze the effect of edge defects [3]. In this model, the 
solar cell is divided into two regions, a central region and a perimeter or 
edge region. The width of the edge region is determined by the low-level 
diffusion length of the minority carriers in the base. A parameter, called 
the. diStCLnce. .cn^^uence, was introduced to characterize the range of a 
defect and to relate the edge region width to the diffusion length. The simple 
one-dlmenslonul low-level diffusion solutions of p/n Junction current-voltage 
characteristics are then obtained for these two regions. The reduction of the 
open-circuit voltage was then c<miputed numerically by connecting the two solar 
cells of these two regions in parallel. 

The idea behind this developed perimeter model and the concept of 
dlstance-of-lnf luence can also be applied to the randomly distributed bulk 
defects of random sizes across the BSF Junction. To Illustrate this defective 
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device model, a partial cell with several unlformaly space defects Is shown In 
Fig. 3.4(a) which has three repetitive square unltcells enclosed by dashed 
boundaries. An enlarged top view of: a unitcell Is shown In Fig. 3.4(b). 

Extension of this uniform defect model to randomly distributed defects should 
become obvious after describing anil analyzing the uniformly spaced defect model. 

The unit cell Is divided Into three regions for analysis purpose. They 
showing the lateral Influence of the defects. They are illustrated In Figs. 

3.4 (b), (c) and (d) and described below. 

(1) Region III 

This is the central region of the unit cell which contains the 
defect. Its cross sectional area, A^, Is equal to the defect 
area, A^. 

(2) Region II 

This Is the first *<tdge region' of the unit cell which surrounds 
the central region. The p/n junction In this region will 
experience the pres<mce of the defective BSF low-high junction 
In region I. It has an area of A^. Its width, W^, Is determined 
by the distance of jLnfluence or the range of the defect, L^, and 
Is given by (Fig. 3.4(c)] 

W 2 - /lJ - (1) 

where Is the thickness of the base layer. 

(3) Region I 

This Is the second 'edge region' of the subcell which surrounds 

the first edge region, region II. The electrical characteristics 

of the p/n junction In this region will not be affected by the 

defect In region III since It is outside of the range or the 

distance of Influence of the defect. Its area is denoted by A^. 

This region does not exist If the defect density, Is high or 
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Fig. 4 (a) Top view, (b) expanded unit cell top view, (c) expanded 

cross sectional view, and (d) developed perimenter cross 
sectional view of a defective back-surface-field solar cell 



















the dl8tance*o£-influence is large# The condition for the 


disappearance of this region or is when 


Lj > 

(2) 

or 



(3) 


This three-region device model for a bulk defect differs from the 
two-region device model for the edge defect in having a finite region III 
which contains the bulk defect. In the two- region device model for the edge 
defect, region III has zero area since the edge defect is on the perimeter 
surface of the cell and has zero projection onto the plane of the p/n junction. 

Although analytical three-dimensional series solutions can be obtained 
in the three regions for the low-level dark diffusion current-voltage characteristics, 
simpler solutions will be worked out using the one-dimensional solution and 
developed perimeter model. The approximation in such a one-dimensional model 
is probably not much worse than the analytical three-dimensional low-level 
solution since in a real cell the defects have random size and are randomly 
distributed while the analytical three-dimensional solution needs to define the 
shapes of the subcell and defect. The one-dimensional solution, to be given 
below, is physically transparant and simple to grasp, and allow us to study 
the effect of defect density and size on the cell performance reduction using 
families of curves which can be computed with very little computing time. 
Three-dimensional analytical solution would give much more complex solutions 
which would not be as transparant and require much more computing time to 
generate the same families of curves to estimate the numerical reductions of 
cell performance due to the presence of the defects. 
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Th* CToss^saetional view of the three reflone of the unit cell le ehown 
In Fig. 3.4(c). The developed aodel le shoim In Fig. 3.4(d) where regions 
I and III have constant base thicknesses equal to the actual base layer thickness . 
X^. Region II has a varying base thickness » Increasing from X^ to L^. 

The distance of Influence or the range of the defect, L^, Is assumed to be 
proportional to the minority carrier diffusion length In the base, Lj^. In the 
numerical calculations to be given later. It will be shown that the nunerlcal 
results on the open-circuit voltage reduction will have reached an asymptotic 
and constant value when Is taken to be greater than about 1. 
of the numerical results obtained In this paper are for L^«2Lg. 

The one-dlmenslonal dark diffusion current solution for a p/n Junction 

at low level with a back-surface field was first applied to BSF solar cells by 

Godlewskl, Baraona and Brandhorst [4]. Tt^ils solution is applied to all three 

regions of the defective unit cell. For the defect-free region I, a zero 

Interface recombination velocity is assumed at the base-side of the BSF low-high 

Junction. For the defective Region III, the defect is characterized by a finite 

Interface surface recombination velocity, again on the base-side edge of the 

space-charge layer of the BSF low-high Junction. For the worst defect, which 

corresponds to a defect that completely short-circuits the BSF Junction. It 

will be shown in section 4 that when the normalized interface recombination 

velocity, S_/v. (where v ■D_/L_ Is the diffusion velocity of the minority carrier 
U y J> IS 

In the base layer), is greater than about 10, the open-circuit voltage degradation 
will have essentially reached its highest value and the defect acts essentially 
like a short circuit or the defective Region III Is essentially a BSO (back- 
surf ace-ohmlc) cell. 

In region II, the one-dlmenslonal dark diffusion current solution Is also 
employed. The varying base layer thickness In Region II, In Fig. 3.4(c) and (d). 


SLg. Thus, all 
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is C«lian ear« of by an avaraga proeadura which waa alao aaployad for tha Itaglon 
II of tha adga dafact [3]. In this avaraga* tha ona-dlmanalonal diffusion 
currant danslty at aach basa layar thlcknass* X* la aunoad ovar tha antira 
Raglon II* froB X»X^ to X*L^ along tha aurfaca from y*0 to 
datailad algabra is glvan in the appandix. 

Tha dark currant-voltage charactariatica of tha thraa raglons ara 
suamariaad below. In Regions I and III* the current densities are constant 
across the plane of the Jimction ao that the current in each region la Just the 
current density times the area. In Region II* the varying current density is 
integrated over its area to give the total current whose detailed steps are 


given in the appendix. The dark currents are 

^1 " - 1] (4) 

I3 - JjA3[exp(qV/kT) - 1] (6) 

The dark current flowing in an Ideal cell without any defects is given by 

Iq “ JQAQ[exp(qV/kT) - 1] (7) 


The current densities and the areas of the three regions as well as in the 
defectless cell are given by 

Jq - Ji - (qDj,Pg/L3)tanh(Xj/L3) (8) 

pxp I (Lj/Lg) tanh (Xg/L^) ] + (SpLg/Dg) sinh (L j/Lg) 


J2 - 
*^3 • 


jexp((Xg/Lg)tanh(Xg/Lg)] + (SgLg/Dg)slnh(Xg/Lg)J 
;inh(Xg/Lg) + (SgLg/Dg)cosh(Xg/Lg) 


lcosh(Xg/Lg) + (SgLg/Dg)slnh(Xg/Lg)J 

Aq » Aj^ + A^ + Aj (11) 

A2 - + /Aj7 (12) 

Dg . 1 /Aq (13) 

Dg is the areal dexisity of the defect since Ag is the unit cell area and 
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(9) 


( 10 ) 



ch« unit €«11 Mlnetnd no Uint It contlMia only ono dnfoet no Ulaotrncod in 

th« unifoni nodol shoun In Figs. 3.4(a) sad 3.4(b). Tbs otiMr psraastsrs In 

I 

(8) to (13) have tbs convantlonsl mssnlng. q Is tbs asgnltuds of tbs slsctron 
cbsrgs. P^, and ars rsspsctlvsly tbs dlffuslvlty, squlllbrluB 

t 

concsntratlon* diffusion langtb and llfstlas of minority carrlsrs la tbs quasi* 
nsutral bass laysr undsr low*lsvsl condition. la tbs rangs or dlstaacs*of« 

. Is tbs tblcknsss 

of tbs quasl*nsutral bass laysr. Is tbs sffsctlvs Intsrfacs rsconblnatlon 
vsloclty of a dsfsctlvs low-blgh Junction. Is ussd to rsprsssnt dsfscts 
which ars not completely electric short circuits across the BSF lowhigb 
junction and defects which are due to a hlgh*layer (p4> layer of the p/p^ 
low*hlgh Junction or n^ layer of the n/nf low-high Junction) of very high 
recombination rate or very low minority carrier lifetime. It can also 
represent partial penetration into the high-layer by the back ohmic contact 
metal or alloy as well as diffusion pipes. 

There are three independent parameters which characterize the property 
of the defect: (i) the defect area, ^"^ 2 * defect density, D^*1 /Aq, 

aiul (ill) tv-, surface recombination vsloclty at the defect, The area of 
Region II which experiences the presence of the defect is related to the 
defect area by (12). The area of Region I which Is outside of the range or 
the influence of the defect is related to the defect area and density by (11). 
Vhen a cell is very thin such that its thickness la less that given by (3) , 


Influence of the defect and Is taken as 2L^ In this paper 




p*rf«cc cell b« tlMo at op«n-clrcult or coro load currant, wa hava 

- Ij ♦ I2 + I3 - (J3A^+J2A2‘M3A3){axp(q(V^‘«VQj,)/kTj - 1) (14) 

for tha dafactlve cell and 

*l ■ *0 ■ 

for Che perfect cell. Normally, ^Qg+^V^>>(kT/q) so that we can drop 1 

compared with «xp(q(VQg+aV^)/kTj or exp[q(V^/kT)J In (14) and (15). Then, 

the reduction of is independent of the photocurrent, I^, or the open- 

circuit voltage which can be eliminated between (14) and (IS). This gives 

- - (kT/q)log^t(A3/AQ) + (J2A2+J3A3 >/JqAqJ (16) 

Numerical solutions and families of curves have been computed from the 

general solution given by (16) using (8) to (13) for the current densities and 

the areas. A little algebra of (16) will show chat all the length parameters 

can be normalized to the diffusion length, L^. Thus, the normalized variables 
2 2 

are where Ajj*A3, and DpL^ where D^^I/Aq. The defect surface 

recombination velocity, S^, can also be nocmalized to a minority carrier 

diffusion velocity, Thus, the reductions of Che open-circuit voltage 

are graphed as a function of these normalized parameters. An lnqjortant general 

result is that the reduction of the open-circuit voltage Is independent of the 

area of the cell but depends on the areal density of the defect, D^. Although 

our analysis was made for a periodic unit cell of the uniform defect model, thla 

general result Is valid for randomly distributed defects. In that case, 0^^ 

would be the average density of the defect over the entire cell which may 

vary from one subcell to the next subcell. 

These numerical solutions, plotted as a function of thickness of the 

quasi-neutral base layer (or roughly the cell thickness since the n4> and p4> 

layers of the diffused nt-/p and p/p4* junctions are less than one micron thick), 

2 2 

with the three defect parwneters, and D^, as constants are given 

in the next two sections with their physical interpretations. 
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3. EFFECT OF SHORT-CIRCUIT BULK DEFECT ON OPEN-CIRCUIT VOLTAGE 

We shall first investigate the worst ^.ulk defect which Is the defect that 
completely short circuits the BSF low-high junction over the defective area, 

In this case, and the current densities in Regions II and III given by (9) 

and (10) simplify to 


sinh(L 




3.1 Dependence on the Dlstance-of-Influence 


Our model relies on the parameter L^ which is the range or the dlstance- 
of Influence of the defect. It was assumed In the model that outside of this 
range or dlstance-of-lnfluence, the defect will not affect or Increase the 
dark recombination current and hence nor the open-circuit voltage. Thus, a 
test of the goodness of this model would be to determine how much the results 

will change with different choices of L^. For this test, we assumed a defect area 

2 2-2 
of 100 pm , a diffusion length of 100 pm so that A^/l^ * 10 . We also 

2 

assime a normalized defect density of D L_ » 0.1. The fractional defective 

U D 

-3 

area In this case Is ApD^*10 or 0.1% of the cell area Is defective. 

The reduction of the open-circuit voltage, -AV^, computed from (16), 
at three cell thicknesses, X^/Lg*0.01, 0.1 and 1.0, are shown In Fig. 3.5. It 
Is evident that for rhin cells (X^/Lg^O.Ol and 0.1) the results do not change 
much If the range or dlstance-of-lnfluence Is assumed to be 1.5 times the 
diffusion length or greater. For the thicker cell (X^/L2=1.0), the reduction 
of the open-circuit voltage by the defective area is small (only about 10 mV), 
but In order to get this result, the distance of Influence must take a value 
greater than twice the diffusion length. Thus, In the remaining calculations. 


L^-2L^ Is selected. 
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Fig. 5 The sensitivity of the open-circuit-voltage reduction on 
the choice of the distance of influence. 
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3.2 Dependence on the Defect Area 

The dependence of the open-circuit-voltage reduction on the defect area, 

Aq, or the area of Region III shovm In Fig. 3.4(b) Is Illustrated In Fig. 3*6. 

Here, the reduction of the open-circuit voltage Is given as a function of the 

2 

normalized defect area, for three cell thicknesses, X^/L^-O.Ol, 0.1 and 

2 

1.0. The constant defect density is DjjLg“0.1. We also used and 

The results in Fig. 3.6 show that the open-circuit reduction, - ^^qq* 

reaches constant values when the defect area is smaller than about 0.01 times 

2 

the diffusion length squared, L^. This asymptotic behavior can be understood 

from the geometry of the device model. From Fig. 3.4(b), It Is evident that 

when or the defect area A^ becomes small. Its contribution to the base 

recombination current becomes small. The Increase of the base recombination 

current comes entirely from Region II when A^-K). The relative contribution 

from the two regions, II and III, Is slightly smaller than their area ratio, k^/ky 

since the p/n Junction In Region II is further away from the defect than in 

2 

Region III. When A^ becomes small compared with L^, A^ approaches a constant. 

From (12) , when i/A^ fs about 10 times smaller than /Lj-X^= 2L^/l-(X^/2Lg)^, 

then A 2 is essentially a constant given by 4(Lj-X^). 

This insensitivity to the defect area is an important result since in 

practice the defect area is probably small compared with the diffusion length 

squared in hlgh-efficlency cells where the diffusion length Is substantially 

2 

greater than 10 urn. Thus, In practice, is less than about 0.01. 

The effect of defect area can also be examined from the opposite direction. 


i \ 


Fig. 3.6 shows that open-circuit degradation, - begins to Increase when 

2 

the defect area Is greater than about 0.1L-. This corresponds to a fractional 

defective area of DjjAj^®Oelx0.1*0eOl in this example or 1% of the area is defective. 

2 

This insensitivity is also the basis of taking A^/L^-OaOl in the analysis of 
the dependence on given in the preceding section, 3.1. 
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DEFECT AREA Aq/Lb 

Fig. 6 The reduction of the open-circuit-voltage of a BSF solar 
cell as a function of defect area with cell thickness as 
the constant parameter. 









3.3 Dependence on the Defect Density 


In the preceding section. It was shown that Increased recombination or 
reduction of open->clrcult voltage came mainly from region II, the peripheral 
region that surrounds the defective area. Region III, since in practice, the 
size of the defective area Is much smaller than the diffusion length. This 
Illustrated the lnq>ortance of the peripheral area surrounding the defect. 

As the defect density Increases, the nund>er of defect Increases for a 
given total cell area. Thus, the total surrounding area (the sum of the areas 
of each Region II) also Increases since It Is proportional to the number of 

defect. The dependence of the open-clrcult-voltage reduction, - AV^, on the 

oc 

defect density Is expected to be large. This Is Illustrated In two ways: - AV„ 

OC 

2 2 

Is plotted as a function of defect density, D^L^, with (1) defect area, 

kept constant or (11) fractional defective area, D^A^, kept constant. These 

are shown In Figs. 3.7(a) and (b). The large dependence and. Increase of the open- 

clrcult-voltage reduction with Increasing defect density Is evident In Fig. 3.7(a) 

2 

when the defect area, A^/L^, is kept at constant values. When the fractional 
defective area reaches unity, reaches the value of a BSO cell which has no 
BSF low-high junction and - AV^^ also reaches Its highest value which is the 
difference between the BSF and BSO cells. For the example shown In Fig. 3.7, 
this asymptotic value for - AV^^ is 235.8 mV. 

The general trend of rapid Increase of - AV^^ with Increasing defect 
density Is also evident In Fig. 3.7(b) when the fractional defective area Is 
kept at constant values. The Increase Is somewhat slower than when the defect 
area Is kept constant since as D^ Increases, A^ decreases In order that Is 

constant and there Is a decrease of - AV^^ at large A^. 

The Important practical Implication of this result is that the defect 
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Fig. 7 The reduction of the open-circuit-vol tage of BSF so]'.r cells 
as a function of defect density with defect area as the 
constant parameter (top figure) and fractional defective area 
as the constant parameter (lower figure). 







density should be kept as low as possible In order to obtain the maximum 
benefits from the BSF junction. For two cells each having IX defective area, 
the one with the smaller defect density and larger Individual defect area will 
give higher open-circuit voltage and efficiency. Thus, If a choice must be 
made in a practical situation between defect density and Individual or average 
defect area, better open-circuit voltage and efficiency Is obtained by keeping 
the defect density low and defect area large. A cell with one defect will 
have a hi^er open-circuit voltage and efficiency than a cell with two defects 
of half area. The physical reason is that the one with two defects will have 
two peripheral areas, one surrounding each defect, idilch would give nearly twice 
as much recombination current as the one peripheral area In a cell with one 
defect. 


3.4 Dependence on the Cell Thickness 

The results of the analyses in the preceding three sections, 3.1, 3.2, and 

3.3, enable us to select a range of values for the defect parameters, area and 

density, so that they fall Into practical ranges. The dependence of the open- 

clrcult-voltage reduction on the cell thickness will be Investigated using 

-4 2 

these ranges of values. The defect area Is set at 10 since It was shown 

2 

In Fig. 3.6 that - AV^ Is Insensitive to the defect area when Is smaller 

than about 0.01 and since In practical high-efflclency cells, Is usually 

greater than 30 or 100 pm and the defect area Is not likely to be greater 

2 2 
than 100 pm or 10 pm on a side. Thus, in practice, A^/Lg is likely to be 

2 •A ^2 

less than 0.01. The choice of Ap/Lg«10 ^ Is made since below 10 , - 

Is nearly constant and approaches its asyiq>totlc value at 10 . 

Sg ■ • is assumed and the dependence on this parameter Is analyzed In 
a later section. L^*2Lg since the result becomes Independent of when It Is 
greater than about 1.5 as discussed In Section 3.1. 
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The defect density Is the constant parameter. The range of the normalized 

defect density, Is from about lO"^ to 10“^ for high efficiency cells which 

is estimated as follows. For a high efficiency cell, is between 100 ym and 

1000 ym. We let Lg-100 um. is equal to the number of defect divided by the 

2 

cell area, D ■N-/A-. A production cell has an area of about 50 cm (3 inch 
D D C 

2 

diameter) or greater* We let A *100 cm * For a good cell fabrication process » 
the number of defect is probably considerably less than 1000 and let to 1000. 

Thus, for a hl'jh efficiency cell, the normalized defect density should be in 
the range of -(Njj/A^)Lg- (1/100) (100xl0"^)^-10"^ to (1000/100) (100xl0“^) ^-0.1. 

Thus, for a low defect density cell, the defect density range is approximately 
DpLj“10*"^ to lO”^ while for highly defective cells, DpLg > 10 

The range of the normalized cell thickness, Xg/Lg, covers from about 
10~^ or less to about 10. This range covers a diffusion length of 100 ym 
and cell thickness of 1 to 500 ym. 

From the above estimates of the ranges of the parameters, the reduction of 
the open-circuit voltage as a function of cell thickness is computed from (16). 

_3 

The results for the entire defect density range of 10 to » are shown in the 

upper part of Fig. 3.8. The expanded graph, for low-defect-density hlgh- 

2-4-1 .f 

efficiency cells covering D^Lg^lO to 10 , is given in the lower part of 

Fig. 3.8 with a larger number of the density parameter per decade. 

2 

When the defect density is infinite (curve DgLg«“), the reduction of the 
open-circuit voltage reaches the asymptotic value which is the difference between 
a perfect BSF cell and a BSO cell which has no BSF junction at all. The upper 
part of Fig. 3.8 shows that the asymptotic dependence on cell thickness is 
approached when defect density, Is greater than about 10 to 100. For this 

limiting case, (16) simplifies to 

- AVqj, - 2(kT/q)log^(ctnh(Xg/Lg)] (19) 

which gives the linear dependence on the semi-log plot in Fig. 3.8 for thin cells. 
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A nunerlcal example la given next to Illustrate the possible Importance 
of defective BSF Junction on the open-circuit voltage. For this example, we 

use the parameters of the p+/n/n+ cell which has the highest reported V*_ of 

oc 

about 620 mV [2]. The reported values are t^^AOO ms, 1^"775 Mm, X^*314 pm, 

Xg/Lg«0.4. Suppose that there Is one short-circuiting defect per cm^, then 

2-3 2 
>6.0x10 . Let us suppose that the normalized defect area, Is 

_2 

less than 10 so that the result Is Independent of Its area. This 

corresponds to a defect size of less than 80 pm on an edge, for example, the 

2 

defect area can be taken as 1 pm . Then, > - 5.49 mV from (9) or using 

2 

Fig. 3. 8. If the defect density Increases to 2 per cm , then AV^^ ■ - 10.0 mV. 

These two numerical examples show that reduction of the open-circuit voltage 

In a hlglt-efflclency cell due to defects across the BSF Junction may be quite 

significant even for fairly thick cells. 

V reduction Is substantially larger If the cell Is thlimer. For 

example, maximum efficiency of the above cell would be reached when the cell 

thickness Is reduced to about 80 pm as Indicated In Fig. 1 or X_/Ln * 0.1. 

B B 

2 

Using the above numerical example and supposing that Djj“l defect/cm , then 

AVq^ “ ~ 25.697 mV while for 2 defect/cm^ we have AV^^ “- 38.26 mV. These 

2 

two numerical results show that the presence of 1 or 2 defects In 1 cm of cell 
area would reduce the expected improvement of V^^ of approximately 60 mV by as 
much as a factor of 2 to only 30 mV when the cell thickness Is reduced from 
the conventional 300 pm to an optimum 70 pm. 

3.5 Dependence on Diffusion Length 

In the preceedlng analyses and numerical examples, the minority carrier 
diffusion length In the base layer, L„, Is the normalization parameter for 
length and area. Dependence of V^^ on defect density, defect area and cell 
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thickness wss anslyzed by keeping the diffusion length a constant. In sane 
practical situations, such as debugging a fabrication process, evaluating 
radiation hardness and damage and others, It may be Important to know how 
much the cell performance will be changed by the presence of defective BSF 
Junctions idien the minority carrier diffusion length or lifetime Is changed. 

Instead of presenting yet another family of curves of - vs Lj^, we 

make use of Fig. 3.8 to Illustrate the dependence of on L^. From the 

two parts of this figure. It Is evident that If we keep the cell thickness, X., 

and defect density, D^, constant and Increase L^, then will decrease 

2 

_ niir II 1 

OC 


which would Increase the reduction of but will Increase also idilch will 


further Increase the reduction of or - 

For exai^le, If the diffusion length Is Increased by a factor of 10 from 

2 

L-/X >1 to 10 and O.!.- is correspondingly Increased by a factor of 100, say from 
B B D B 

-1 

10 to 10 , than would be Increased from 0.147 mV to 78.349 mV. This 

would substantially nullify the expected Improvement of when Is Increased 
or the minority carrier lifetime is improved by some processing or annealing 
procedures. 

On the other hand, If a cell has some defects across the BSF junction, 

loss of lifetime or diffusion length due to radiation or other reasons may not 

decrease the open-circuit voltage as much as a cell which has no defects across 

the BSF Junction. This Is for the reason that Is already lowered In the 

defective BSF cell and decrease of L_ would make the defective BSF areas less 

B 

effective In reaching the front Junction to Increase its recombination current. 
This Interesting result Indicates that one could Increase the radiation 
resistance or make V_ less sensitive to radiation Induced recombination center 
over a wider range of radiation dose when some localized defective areas are 
Introduced Into the BSF Junction which would reduce Initially but also 
stabllze It with respect to the change of L^. 
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4. EFFECT OF BULK DEFECT Wll^ IMTERMEDIATE RECOMBINATION VELOCITY ON 

OPEN-CIRCUIT VOLTAGE 

In practlC8( the defects across the BSF junction may not be completely 
a.'iort-clrcultlng the BSF junction electrically. In addition* some defective 
BSF junction areas may contain a high concentration of dufects or recombination 
centers tfhlch are distributed In the entire layer of the hlgh-Iayer (n^ layer 
of a n/n4- BSF junction or p4 layer of a p/p*** BSF junction). These and other 
defective BSF junction areas can be represented by an effective surface 
or Interface recombination velocity across a defective area, Sj^. This was 
Included In the current densities for Regions II and III given by (9) and (10) 

In section 2. S^^ Is the effective surface recombination velocity of minority 

carrier at the boundary plane on the base side of the space-charge layer of the 
low-high Junction. It Includes all the recombination In the space-charge layer 
of the low-high junction as well . the recombination In the high-layer and on 
the contact surface of the back metal-semiconductor contact. 

It Is evident from (9) and (10) that the proper normalization for Is 
Dg/Lj"/Dj/Xg which Is the minority carrier diffusion or recombination velocity 
In the base layer. Physically, we can consider the two rec(»Dblnatlon processes 
In the two regions, base and BSF regions, as two serial processes. If base 
recombination Is large, then It dominates. If Xs large, we have two possible 
situations: thin and thick base. For thin base, then Sj^ dominates 

but the recombination current Is determined by base thickness. For thick base, 
base recombination dominates and will have little effect. 

The dependence of the reduction of the open-circuit voltage on the BSF 

defect surface recmnblnatlon velocity Is Illustrated by a family of curves 

/ ^ 

shown In Fig. 3.9. In this figure, Ajj/Lj«10 , DpLj*0.1, Lj*2Lj and 

Is the constant parameter with a range of. 0.01 to 100. The reduction of the 
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-AV^-(mV) 



10 '^ ^ 10 "' ® 10 

THICKNESS Xb/Lb 


g. 9 The reduction of the open-circuit voltage of BSF solar 
cells as a function of cell thickness with the defect 
surface recombination velocity as the constant parameter. 






op«o»clrcule voltage, - ^^OC* easentially reached lea limiting value tdien 
the defective aurface recombination velocity la greater than ten tlmea the base 
diffusion velocity or > lOD^/I^. 

In a praurlcal high efficiency cell, Dg*10 cm^/sec and I^^lOO to 1000 ua 

or 0.01 to 0*1 cm so that * 1-000 to 100 cm/ sec. Thus, the limiting or 

3 4 

highest *> Is reached when Is greater than 10 to 10 cm/ sec. This 

range of Is readily attained and exceeded when physical defects are 

present across the BSF Junction. Smaller values are given by defects which are 

more spread out Into the p4> or n4- layer of the low-high junction. 

The main point shown in Fig. 3.9 Is that reducing the defect surface 

recombination velocity will not improve or decrease reduction since 

2 

it is necessary to reduce to less than 100 cm /sec to have any significant 
gain. It Is more effective to reduce the defect density or the number of 
defect than to reduce the defect surface recombination velocity since the base 
diffusion-recombination velocity of a high efficiency solar cell is rather low, 
100 cm/sec to 1000 cm/sec and since In practice. It Is difficult to achieve 
surface recombination velocity of values less than about 100 cm/ sec for 
Interfaces and bare or oxidized surfaces and probably more difficult to achieve 
values equal to or less than this for a defect surface. 

The results of bulk defect can also be compared with that of edge defect 
treated in another paper (3]. At a given Sj^/CD^/L^) such as 1.0, bulk defect 
reduces more than edge defect. Consider a small cell to emphasize edge 
defect. For the same reduction of V^, of edge defect must be more than ten 
times that of a bulk defect, again because bulk defects have much larger 
peripheral areas (Region II) than edge defects. 






5. SUMMARY AMD CONCLUSION 

Hlgh-'efflclency solar cells contain a low-high hack-surface-fleld (BSF) 
Junction (n/nf or p/p4-) to provide a large potential barrier «^lch prevents 
the photogenerated and front-junction Injected minority carriers from reaching 
the hlgh-recomblnatlon-rate ohmic contact on the back surface. Such an 
isolation greatly reduces the recombination current and substantially Increases 
the open-circuit voltage and energy-conversion efficiency. The low-high junction 
or the back-surface field layer is very thin (about 1 ym) and can easily be 
shunted electrically by material and manufacturing defects introduced during 
crystal growth and cell fabrication processes. This paper gives a mathematical 
analysis of the effect of defects* across the back-surface-fleld junction and 
located away from the cell edge* on the solar cell performance. Analysis of 
the edge defect is given In another paper [3]. 

In this paper* several families of design curves are presented which show 
how much the open-circuit voltage of a back-surface-fleld solar cell Is reduced 
as a function of the cell thickness* defect area* defect density and Interface 
or surface recombination velocity at the defective area. 

The defects across the BSF junction are randomly distributed and have 
random shapes and sizes which defy exact analytical or numerical solution. 

The random defects are modeled by an uniform device model which Is characterized 
by three defect parameters: the average defect area* the average defect density 
(per unit area) and the effective surface recombination velocity of the defective 
area. A dtveZoped pZAAmttdA. device model Is employed in order that the 
approximate one-dimensional analytical diffusion solutions of p/n junction with 
BSF can be applied to a defective BSF solar cell. In this model* a defective 
cell with many defects Is divided Into many one-defect unit cells. A one-defect 
unit cell is further divided into three regions (Fig. 3.4): a central region* 
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Region III, which Is the projection of the defect onto the front p/n junction 

plane, and two concentric surrounding or peripheral regions. Region II and I. 

The one-dimensional solutions are applied to Regions I and III. The cross- 

section of the middle region. Region II, Is de.vziopejd 30 that the one- 

dlmenslonal solution can also be applied to this region of varying thickness. 

Its width Is determined by a distance parameter, L^, which Is the range or 

the distance of Influence of the localized defect. It Is the nuiximum distance 

at which the p/n junction Injection-recombination current can still be Influenced 

by a localized defective BSF junction area. This distance of influence Is 

shown to be about two diffusion lengths of the minority carriers in the base 

region (Fig. 3.5). The main results of this paper are summarized below. 

The reduction of the open-circuit voltage is nearly Independent of the 

2 

defect area when it is smaller than about O.IL^ or when its dimension Is less 

than about 30% of the minority carrier diffusion length in the base (Fig. 3.6). 

This insensitivity on defect area, when the area is small, arises from the fact 

that a defect will increase recombination and reduce open circuit voltage not 

only on the front junction area which is the projection of the defective BSF 

area (Region I) but also the junction area surrounding the projection (Region II). 

As the defect area decreases, the peripheral recombination area (Region II) will 

stay relatively constant, giving rise to a defect-area-independent ^OC reduction. 

In practice, there should be little defect-area dependence since in good cells, 

2 2 

dei.ects are mainly small size and not likely to be as large as O.lLg^O. 1x100 

2 

^1000 urn . This also means that small-area defects will be as damaging as large- 
area defects. 

The main influence of the open-circuit voltage reduction comes from the 
density of the defect (Fig. 3.7). The reduction of Increases with 
increasing defect density since there are more peripheral recombination areas, 
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one per each defect » when the number of defect In a cell Increases. For 

example, if the fractional defective area Is kept constant (such as the 0.01 

or IZ curve of the bottom graph of Fig. 3.7) and If a single defect of 100 um^ 

Is divided Into 100 separate defects of 1.0 each, the reduction of the 

open-circuit voltage will Increase substantially because the total peripheral 

recombination area (sum of Region II) will have Increased 100 folds. 

This general trend will persist even when three-dimensional effects and high- 

level condition are taken Into account as the defect area becomes small and 

the defect current density Increases. 

Because of the importance of the perimeter region whose size Is given 

by the range or dlstance-of-lnfluence, the presence of even one defect of 

small area across the BSF junction will cause significant reduction of open- 

circuit voltage of a hlgh-efflclency BSF solar cell. A numerical example 

2 

showed that If there Is 1 defect per cm of cell area In a 320 pm thick BSF 

state-of-the-art hlgh-efficiency cell (L_=775 ym, V *620 mV and EFF=17%), 

D OC 

one would expect a reduction of open-circuit voltage of 5*5 mV. If there 

2 

are 2 defects per cm , the reduction of would increase to 10 mV. 

The reduction of due to defects across the BSF junction increases 
when the cell becomes thinner. This was the initial consideration that 
motivated this analysis since high-efficiency BSF cells reach their peak 
efficiencies when the cell thickness is decreased to about O.IL^. The larger 
reduction of in thinner cells is again due to the lateral Influence of 
the defects which increases recombination in peripheral region that surrounds 
the defect (Region II in Fig* 3.4) since this regl^^n Increases in width when the 
cell gets thinner. As a numerical example, if the state-of-the-art hlgh- 
efficiency cell just described has a thickness of 80 ym instead of 320 pm, 

2 

then the reduction of is increased from 5.5 mV to 25.7 mV for 1 defect/cm 
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2 

and from 10 mV to 38.3 mV for 2 defects/cm . 

Families of design curves on the thickness dependences of V^^ are given 
In Fig. 3.8 with the defect density as the parameter. 

Defective areas across the BSF junction may not be complete electric 
short circuits but shunts with non-zero resistance. This was represented by 
an Interface or surface recombination velocity to characterize the defective 
area. The effect of the surface recombination velocity on V^^ reduction was 
analyzed which showed that V^^ Is Improved only when It Is substantially smaller 
than ten times the diffusion velocity of the minority carriers In the base, 

Dg/Lg (Fig. 3.9). For hlgh-efficlency cells, the diffusion velocity Is not 
more than about 100 to 1000 cm/ sec. Thus, to have any Improvement, the defect 
surface recombination velocity must be made substantially smaller than 100 
to 1000 cm/sec. In practice, a defective region rarely has recombination rate 
and velocity less than about 1000 cm/sec. Thus, a defective region acts 
essentially as an electrical short circuit across the BSF junction, otherwise 
it is a good area. There are few inbetween. 

From the results of this analysis, it is evident that material and 
fabrication defects across the BSF junction of high-efficiency cells can be 
the major deterrent to high open-circuit voltage and high efficiency even if 
only a few or even one defect is present across the back-surface-field junction. 


6 . 


APPENDIX 


6.1 Analysis for Region II Surrounding the Defect 

The analyses which led to the current density expressions given by 
equations (8), (9) and (10) in the text are given In this appendix. For 
Regions I and III, the current density equations, and J^, are the standard 
solutions of a p/n junction with a BSF junction whose surface recombination 
velocity Is 0 (Region I) and (Region III which has the defect) . These 
solutions are one-dlmenslonal low-level diffusion solution in the quasi- 
neutral base region of a p+/n/n+ or a n+Zp/p-f diode with a constant 
bulk Impurity concentration. The surface recombination velocity, S^, of the 
defect In Region III can be thought of as an effective recombination velocity 
if a refined model Is needed to represent defects which are extended Into the 
nf or p+ layer of the BSF junction or point defects or small defect clusters 
which are distributed In this layer. In such cases, Sq is given by [4] 


[8inh(X./L.) + (S.L /D.)cosh(X,/Lj] 

S_ - (D./L.)(P./P_) ^ 

[cosh(X^ZL^) + (S^L^/D^)sinh(XyL^)J 


(A.l) 


where D_j_, and are the diffuslvlty, diffusion length and concentration of 
the minority carriers in the heavily doped n+ or p+ layer. X^ is the thickness 
of the n-f or p+ layer and is the surface recombination velocity of the 
ohmic contact on the back surface of the semiconductor. P^ is the minority 
carrier concentration in the quasi-neutral base region. Because of the very 
large ratio of Pb/P. since the BSF layer is much more heavily doped than the 
base layer, S^ is quite small due to the reduction by the factor (P^/P^). This 
factor represents the Boltzmann factor of the potential barrier height across 
the iow-high junction. As a numerical example, D_j_*l cm /sec, 1 ym, 
(D7LJ-10^ cm/sec, P «10^°/10^^-10 cm"^ and P -10^°/10^^«10^ cm"^. Then, 

(D. /L.) (P, /P«)*10^ (10/10^)“! cm/ sec. The fraction in (A.l) involving the 

T + *T* B 

hyperbolic functions is of the order of 1 so that S^ cannot be substantially 
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larger than 1 aa/aec. 

If there is substantial recombination In the space charge region of the 
lov>hlgh junction, then the total would be the sum of that given by (A.l) 
and the effective surface recombination rate of the space charge layer, 
which Is approximately given by thickness of the space charge layer divided by 
the lifetime, Sg^-Wg^/tg^-c^^N^Wg^-Cg^N^g where c^^ Is effective capture rate 

3 

of electrons and holes In the space charge layer (cm /sec), is the volume 
density of the recombination centers In the space charge layer, and Is the 

areal density of the recombination center averaged over the entire volume of 
the space charge layer and Wg^ Is the thickness of the space charge layer. 

Sgg Is also not too large, for example, if Wg^-l pm, ys, then, Sg^«100 cm/s. 

Thus, these examples show that these distributed defects In the BSF layer 
and BSF junction can essentially be considered as perfect BSF junctions since 
their effective surface recombination velocity Is less than about 1000 cm/sec 
and they will have rather small effect on the reduction of the open circuit 
voltage. The large reduction of comes from more serious defects which are 
essentially short circuits across the BSF junction. 

The analysis for Region II is aided by the dzvttopzd p£Avn‘ZteA model shown 
In Fig. 3.4(c) where x is the coordinate perpendicular to the p/n junction and 
y is the coordinate along the direction joining the nearest neighbor defects 
In the plane of the p/n junction. X Is the variable thickness of the base of 
Region II and is a function of y, X=X(y). The total current in Region II is then 


Wo 


^ 2 “ 


Jjj.(X,y)dy Zj.j 


(A. 2) 


where Z is the effective length or the circumference of Region II while W» is 
XX ^ 

Its width given by (1) which Is repeated below. 

(A. 3) 
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The one-dlAenslonal low-level dark diffusion current density Is given by 


Jii(x,y) 




[slnh(X/Lg) + (SjjLg/Dg)cosh(X/Lg)] 
(cosh(X/Lj) + (SjjLg/Dg)slnh(X/Ljj)] 


(A.4) 


where D^, and are the dlffuslvlty, diffusion length and lifetime of 

the minority carriers In the quasl-neutral base at low Injection levels. Is 
the excess minority carrier concentration at the base-side edge of the p/n 
junction and Is given by 


Pq » Pg[exp(qV/kT) - 1] (A. 5) 

where P^ Is the minority carrier concentration In the base, q Is the magnitude 
of the electron charge, k Is the Boltzmann constant and T is the cell temperature. 
Sq was defined In (A.l). 

In order to get an explicit analytical solution for the total current 

flowing In region II from (A. 2), a linear approximation Is made to the quadratic 

relations between the variable base thickness, X, and the position y, so that 

(A. 2) can be Integrated analytically using the current density expression given 

2 2 2 

by (A.4). The exact quadratic relations Is given by X = X^ + y as evident from 
the geometry shown In e'lgs. 3.4(b) and (c). The linear approximation Is given by 
X » X « Xg + [(Lj.-Xg)/»^:pxJ] y (A.6) 

so that the end points are correct and correspond to those shown in Fig. 3.4(d); 

X = x(y«0) = Xg 

and 

X = x(y=W2=/L^-X^) = Lj 


The difference between the exact quadratic and the approximate linear relations 
was demonstrated to be small In analyzing the edge defects [3], giving a maximum 
error of less than 14%. Similarly small errors would be expected here for 
bulk defects across the BSF junctions. 
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The analytical solution of using the linear approximation of (A. 6) 
and the current density given by (A. 4) in the Integral (A. 2), la 


‘2 ■ 


cosh(Lj/Lg) + (SpLj/Dg)slnh(Lj/Lj) 


cosh(Xg/Lg) + (SpLg/Dg)slnh(Xjj/Lj) 


(A. 7) 


The asymptotic solution of (A. 7) for approaches the correct limit, 
however, the solution for is incorrect since in this case, there would be 

no reduction and (A. 7) should reduce to the form of a constant or 


constant base layer thickness which are 

Jll(X-Xj.y) ■ q(DjP„/L,) WnhOj/L,) 


and 


I^CSp-O) 


0<“bW *2 


(A. 8) 


(A.9) 


where A^ was given in (12). But, (A. 7) does not reduce to (A.9). 

In view of the approximations already made, we shall remedy this by 
another approximation so that the limiting solutions are correct. This is 
achieved by replacing the two cosh terms in (A. 7) by exponential terms so that 
they give the correct solution o*" (A.9) when Sjj-0. The replacements are: 
cosh(Lj/Lg) replaced by exp( (Lj/Lg)tanh(Xg/Lg) ] and co8h(Xj/L^) replaced by 
exp[(X„/L_)tanh(X_/L_)]. The new approximate solution for the total current 

DO DO 

in Region II which gives the correct asymptotic result when S^«0 is then 


log. 


exp[(Lj/Lg)tanh(Xg/Lg)] + (SpLg/Dg)8lnh(Lj,/Lg) 


expt(Xg/Lg)tanh(Xg/Lg)] + (SpLg/Dg)sinh(X^/Lg) 


(A. 10) 


This can then be simplified to the form given by (5) with J 2 defined by (9), 
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In Figure 3.4 we have illustrated the developed perimeter model of a 
square unit cell which contains a square-sliaped defect. We shall show that 
the results obtained in the text are also applicable to a square unit cell 
which contains a circular-shaped defect if the width of Region II is redefined 
and allowed to overlap into the four adjacent unit cells. The results listed 
in (11) » (12) and (13) for the areas of the three regions of a square defect 
are derived as follows. From Fig. 3. 4(a) and (b), it is evident that 
A 2 +Aj»( 2 W 2 +*^ 2 )^ for a square defect of area Aj and side where 
is the width of Region II. Since the unit cell area, A^, is given by the 
sum of the area of the three regions, A^^A^'fA^'^A^ and the defect density is 
given by D^>1 /Aq, then the areas of Region II and I can be expressed by the 
defect density, D^, and the defect area A^ which will use the symbol A^. ^us, 
and A^{2\H^^^)^-A^ or 

Ai - - Ajjll + 2(W2/»^)]^ - Ajj (A.ll) 

A 2 « AjjU + 2(W2/V^)]^ - Ajj (A.12) 

”2 " 

(A.12) reduces to (12) if W 2 given by (A. 13) Is substituted into (A.12) and noting 
that we wrote A^^A^. The length of Region II, is defined by 

may be equated to (A.12) of (12) to give '^)“A(/L^-X^ + *^)* In the 

final results given by (4) to (13), does not enter explicitly. 

For a circular or round defect in a square unit cell, the results are 
readily obtained in the same way as above and arc given by 

Aj - D“^ - Ap[l +/S’(W2/v^)]^ - Ajj (A. 14) 

^2 " +»^(W2/*^)1^ - Ap 


(A. 15) 



The factoT» 2 , In (A. 11) and (A. 12) for a square defect is replaced by ■ 1.77 

for a circular or round defect shown in (A. 14) and (A. 15). Thus* If we 
redefine the width of Region II, W^, by *^7 a V ?2 ■ O.SbbW^, then the areas 
A^ and A^ for the round defect have the same expressions as those for the 
square defect. 

6.3 Randomly Distributed Defects 

The developed perimeter device model for a defective unit cell shown 
In Fig. 3.4 for a square-shaped defect and a square unit cell can be extended 
to randomly distributed defects with random distribution of defect area. The 
results obtained In the text In terms of the defect density and area are still 
applicable If these are average density and area over the random distributions 
and If there is a sufficient number of defects so that statistical average Is 
meaningful. In practice, there are few defects In hlgh-efflclency silicon 
solar cells manufactured under controlled conditions so that the discrete 
model of one or few defects uniformly distributed which was described in 
Fig. 3. 4 should be a better approximation than a randomly distributed model of 
many defects. 

Two brief and simple illustrations are given to show how the device model 
can be extended to random defects. It should be noted first that the extension 
to random defects Is considerably simplified using the three-region defect model 
since It is developed to use the one-dimensional solutions in all three regions. 
Thus, the only random parameters that need to be averaged statistically are the 
areas of the three regions, A^, A^ and A^. These areas do not Interact with 
the current densities and the minority carrier distributions in these three 
regions due to the one-dlmenslonal approximation which would have required the 
solution of the three-dimensional diffusion equation with random boundary 
conditions. 
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Consider first the simplest case in which the defect area« is 
varied but still maintains a square shape and constant density or uniform 
separation. In this case, a statistical average needs to be carried out 
for the areas A2 and A^ or A^ over the random distribution of A^^ using 
the expressions given by (A. 11) and (A. 12) or (11) and (12) for A^ and A^ 
which are functions of Aj^. The averaged areas, A^, A^ and A^«Ap, are to be used 
to compute the total current in the three regions using (4), (5) and (6), and 
the open-circuit reduction using (16). A^ la given by D**^ -> A2 ** so that 
only A2 needs to be evaluated. From (12) or (A. 12), it Is evident that we 
need to compute the statistical average of «^. If we assume that the size 
or the edge of the defect, Z^/ky is randomly distributed, then » Z, and 


Ajj-^*(Z3)2. 

Consider next a more complicated example In which the defect spacing 

is randomly distributed. Let the defect area be constant and square and let 

their sides be parallel to the x or y axis. Then, the area A2 flven by 

(A. 12) or (12) needs to be modified since Region II from adjacent square unit 

cells may overlap and counced twice. This overlap area must be subtracted out 

of (12) or (A. 12). To get an expression of the random size of the overlap area, 

let us denote the center-to-center spacing between the adjacent unit cells 

by S and S along x and y directions and S and S are randimt variables. Tnen, 
' X y X y 

there will be no overlap when S^> and > 2W2+Z3 where Zj**^^**^. 

Denote this non-overlap separation by Sj^-2W2+Z3, then the overlap area between 
the two adjacent square unit cells is given by A22* non- 
overlap area of Region II for each unit cell is then ^2~ ^12' spacing 

is still smaller such that some of the overlapped area spreads into Reglcii 
III, then it must be subtracted out of A22 or the net non-overlap area of 
Region II is A2-A22+A32 where A32 is the area of Region III which la • * tapped 
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by th« Region II Iron an neighboring unit cell. The area which nuat be used 
to compute the current flowing In Region II using (5) Is now the statistical 
average of ^2~^22^^32 random distribution of and which vary 

from 0 to « or to the dimension of a cell. 

From these two examples, the procedure for extension to more general 
random distributions of defect density and area as well as shape Is evident* 
although the algebra Is somewhat complicated. They also demonstrate the 
simplicity of treatment for random defects from the use of the three-region 
developed perimeter device model for the defect. 
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IV. THICKNESS DEPENDENCES OF SOLAR CELL PERFORMANCE^ 

1. INTRODUCTKm 

Decreasing the semiconductor thickness of Junction solar cells not only 
reduces the semiconductor material cost but also Improves the efficiency If the 
cell has a back-surface field and Is not too thin: - not less than about one 
average optical absorption length ('V' 10 urn in SI) or about one tenth of the base 
carrier diffusion length In hlgh-efflclency cells (>161 AMI in SI). The improve- 
ment In thin cells arises from the smaller minority carrier recombination volume 
when the cell Is thin. Large number of curves on the performance of solar cells 
as a function of thickness have been computed and given by Hovel (1] for cells on 
silicon and other semiconductors. The performance parameters presented Include 
the open-circuit voltaget the short-circuit current density^ fill 

factor or curve factor* FF, and the efficiency* EFFi, for both space (AMO) and 
terrestrial (AMI) solar Illuminations. 

The theory used by Hovel was the analytical one-dimensional low-level p/n 
Junction diode theory* including emitter space-charge-layer recombination current. 
Because of the low-level assumption* some Important features of the dependences 
of these performance parameters on the thickness of the cell are missing in his 
faisllles of curves. Two most obvious ones are: (1) base resistivity loss when the 
cell Is thick (above about 100 urn), when the resistivity is high and when the solar 
Intensity and cell current are high (concentrator photovoltaic energy conversion) 
and (11) emitter and back-surface-field layer recombination when the cell Is very 
thin. These factors are related to high injection level effects which cannot be 
taken into account readily using low-level analytical theory of minority carrier 
diffusion and drift In p/n Junctions such as those employed by Hovel and later 
workers. In addition to these two effects not taken Into account by the low-level 

# This chapter will be published as an article In Solid-State Electronics In 
1982 with the same title. Zt has been accepted for publication. 
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theory, there are some additional hl^-level effects which are dependent on 
specific materi'al properties, such as the variation of the dlffuslvlty and 
lifetime of the minority carriers In the base layer with Injection or Illumination 
level. The purpose of this short note Is to demonstrate the importance of these 
two universal effects on the performance parameters. In particular, the fill 
factor and the efficiency. 

2. DEVICE MODEL 

The one-dimensional device equations were numerically solved using the trans- 
mission line technique [2]. The Input parameters are the dopant Impurity concen- 
tration profile, the variation of the electron and hole mobility with dopant con- 
centration, the concentration and the thermal capture and emission rates of electrons 
and holes at the recombination center, the develce temperature, the AMI solar Irra- 
dlance spectra and the silicon absorption coefficient. Complete transmissions at 
the front and back surfaces of the cell are assumed. A model recombination center 
of zinc Is used although the results are generally valid for other recombination 
centers and for assumed minority carrier lifetimes. The latter Is a procedure 
most frequently used by solar cell theorists. The device structure parameters 
used in cur numerical analysis are Identical to chose used in reference [2]. 


3. ANALYSIS OF RESULTS AND COMPARISONS WITH LOW-LEVEL THEORY 

The results of our calculations for the Intrinsic cell [3] are presented In 

Figs. 1 to 4 which contain open-circuit voltage, short-circuit current fill factor 

and efficiency as a function of cell thickness from 10 to 1000 ym. Both a 

back-surface-field (BSF) n+/p/p+ and a back-surface-ohmic (BSO) n+/p cell 

12 3 

structure are computed. The recombination center density selected, 10 Zn/cm , 
gives a minority carrier diffusion length of about 577 ym. This was chosen to 
simulate current state-of-the-art high-efficiency silicon solar cells whose 
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residual recombination center could be an impurity-vacancy complex. However, 
the minority carrier diffusion length of a production hlgh-efficlency cell 
will be similar to that assumed in this numerical example. Tlie resistivity 
of the p-type silicon used in this example is about 25 ohm-cm or a dopant 
impurity concentration of SxlO'*’^ boron/ cm"*. 

The open-circuit voltage and the short-circuit current are given in 
Figures 1 and 2. These curves can be compared with those given by Hovel 
based on the one-dimensional low-level analytical theory. For ease of comparison. 
Hovel's curves of Jg^, FF and EFF are reproduced here in Figure 5. It 
is evident that our exact numerical solutions for and J have thickness 
dependences very similar to those giver by Hovel, shown here as Fig. 5(a) and 
5(b). The magnitude difference between Hovel's and our [compare 
Figure 1 with Figure 5(a)] comes from the differences in the solar spectra used 
since his AM2 spectrum has a slightly lower intensity than our AMI spectrum, 
given in reference [2]. The resistivity (10 ohm-cm vs our 25 ohm-cm) and 
diffusion length (232 ym vs our 577 ym) differences contributed to the 
differences in The main feature is that the thickness dependences of 

these tvro parameters, and are not very different using the low-level 

(Hovel) and the exact (our) solutions. 

The similarity ceases for the other two performance parameters, the 
fill factor and the efficiency. The differences in the fill factor is 
particularly important in pointing out the series resistance effect and the 
high injection level effect which are the principal causes of low efficiency 
in thin film cells. The fill factor and the efficiency are given in Figures 3 
and 4 which may be compared with the low-level results given in Fig. 5(b) and 
5(c). 
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THICKNESS (|uM) 


Fig. 1 The short-circuit current density at AMI in a BSF silicon 
n+/p/p-l- and a n+/p solar cell as a function of cell 
thickness. The total available photo current is also 
given. These are computed from the exact solution of the 
Shockley Equations. The base diffusion length is 577 ym. 
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Thickness H. p 


Thickness H 


Thickness H 


The short-circuit current, the open-circuit voltage, the 
fill factor, and the efficiency at AMI in silicon BSF and BSO 
solar cells as a function of cell thickness computed f^z.a 
the one-dimensional low-level analytical theory by Hovel [1]. 






Coaparlson of Fig. 5(b) with Figure 3 shows that there are marked 
differences In the fill factor computed from the low-level approximation and 
the exact numerical solution. The three main differences are listed below. 

(1) The BSF cell (nf/p/pf) has a lower fill factor than the BSO (n-t*/p) 

cell when the thickness is greater than about 25 ym or about 0.043L. 

B 

(Lg*577 ym). The low-level approximation predicts just the reverse 
over the entire thickness range: FF In the BSF cell Is always 
higher than the BSO cell. 

(2) The fill factor peaks between 100 and 300 pm from the exact theory 
and drops rapidly when the cell thickness Increases about 400 pm 
or 0.5Lg. The low-level approximation shows that the FF flattens 
out to a constant and ideal diode value when the thickness Is 
greater than 200 pm and maintains this constant even at Infinite 
thickness. 

(3) The fill factor of the BSF cell Increases towards the ideal diode 
value when the thickness is less than about 25 pm while that of 
the BSO cell decreases continually as the thickness decreases. 

The low-level theory does not show a significant reversal of the 
FF in the BSF cell while in the BSO cell it decreases continually 
when the thickness decreases. 

These differences can be accounted for by the majority carriers and high-level 
injection which are neglected in the low-level minority carrier diffusion theory. 
The differences given in (1) is due to the fact that high injection level 
condition sets in in the base of the BSF cell earlier than the BSO cell because 
the photogenerated minority carrier density in the base of the BSF cell is 
higher while it Is lower in the BSO cell because of the high recombination rate 
at the ohmic contact of the BSO cell. This difference cannot be accounted for by 
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refinement of the low-level theory short of Including the high level effect. 

The difference stated in (2) Is due to the large series bulk resistance of 
the majority carriers when the cell Is thick. In the exact numerical solution, 
this bulk resistance Is automatically Included while In the low-level minority 
carrier diffusion theory, It is not Included since the majority carriers are 
excluded In the lotr-level analysis. The low-level theory can be refined to take 
Into account of the bulk series resistance very approximately If an equivalent 
resistance is added in cerles with the Ideal one-dimensional cell and this 
external series resistance is made proportional to the cell thickness. 

The Increase of the fill factor of the BSF cell towards the Ideal diode 
value (above 0.8) In Figure 3 as the cell thickness decreases (Item 3) Is due to 
the Increasing dominance of emitter recombination when the cell Is very thin since 
base recombination volume decreases with cell thickness. In our example, the 
emitter thickness Is fixed at 0.25 ym. The reason for reaching the Ideal diode 
value Is that the emitter is heavily doped and hence almost always at the low- 
level condition of an Idea diode. 

For the BSO cells, the fill factor decreases with decreasing thickness In 
both the exact and approximate theories [Fig. 3 and Fig. 5(b)]. This Is due to 
the decreasing open-circuit voltage as the cell becomes thinner. One can readily 
show that the FF values given In these figures are nearly equal to those predicted 
by the ideal low-level base recombination theory, Ig^“lQexp(qVQ^/kT) . 

There are also differences In the AMI efficiency obtained by the exact 
numerical and the approximate analytical theories, particularly when the cells 
are thick. This is Illustrated by. comparing Figure 4 with the 10 ohm-cm curves In 
Fig. 5(c). For thick cells, the low-level theory predicts a constant efficiency 
while the exact nisnerlcal solution shows that the efficiency decreases continually 
with increasing thickness. This Is due to the series bulk resistance effect idilch 
is not taken Into account In the Ideal analytical theory. 


- 82 - 


4 


REFERENCES 


1. Harold J. Hovel, Solar Cells . Academic Press, 1975, Chapter 5. 

2« C* X* Sah, P« C« H« Chao, C* K* Wan^, R« L* Sah, K« A* Yamalcawa, 

R. Lutwack, "Effect of Zinc Impurity on Silicon Solar-Cell 
Efficiency," IEEE Transaction on Electron Devices, ED-28, 304-312, 
March 1981. 

3. Series bulk resistance Is automatically Included In the exact one- 
dlmenslonal Intrinsic cell model enq>loyed by us. Contact resistance 
Is assumed to be zero. Lateral resistance loss due to front metal 
grid lines Is a two-dimensional effect not accoutable by the one- 
dlmenslonal model used In our cne-dimenslonal model. 

5. ACKNOWLEDGEMENT 

The author would like to thank Fred A. Llndholm for several suggestions 
to Improvement the manuscript and an anonymous reviewer for Solid-State 
Electronics whose detailed comments were taken Into account in the 
revised manucrlpt given here. 


V. 


THERMAL RECOMBINATION AND GENERATION PROPERTIES OP ELECTRONS 
AND HOLES AT THE DOUBLE DONOR TITANIUM CENTER IN SILICON 


1. INTRODUCTION 

To achieve high energy-conversion efficiency at a reasonable price In 
silicon solar cells, the densities of the residual metallic Impurities In 
the starting silicon material must be reduced at a refinement cost which is 
consistent with the price-performance economics requirement of silicon solar 
cells In Cisapetltlon with other energy sources. High recombinHtion rates 
of electror<s and holes at the residual Impurity centers In silicon will severely 
limit the highest efficiency that is achievable in a solar cell containing 
these residual impurities. The starting metallurgical grade silicon material 
contains many metallic Impurities at high concentrations. These iiq>urlties 
must be removed by chemical refinement processes before single silicon crystals 
can be grown. The various crystal growth processes can also remove a signifi- 
cant fraction of the remaining impurities but It is more economical to remove 
most of the residual metallic Impurities by a prior chemical refinement process 
which Is the proven and traditional purification sequence In the production of 
integrated circuit grade silicon single crystals. 

The Importance of achieving low residual metallic Impurity concentration 
In order to reach the highest energy conversion efficiency in silicon solar 
cells Is only recognized recently [1]. In fact, the purity requirement on 
solar grade silicon for high efficiency performance (>17% AMI) Is even more 
severe than the highest quality silicon single crystal required for the current 
and future hlgh-density very-large-scale high-performance Integrated circuits. 
Thus , It Is Important to determine accurately what Is the aosolute concentration 
of an Impurity that can be tolerated in a hlghest-eff iclency silicon solar cell 
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design when the efficiency of the cell Is specified by the economics requirement 

of a pratlcular appllcatloni such as terrestrial or apace power generation. 

In order to provide this information, the electronic properties of the residual 

Impurity in a single crystal silicon must be accurately known so that a first- 

principle computer calculation can be carried out for a particular solar cell 

design to determine the relationship between the cell performance (such as 

short-circuit current density, the open-circuit voltage, fill factor, 

FF, and maximum energy conversion efficiency, EFF) and the concentrations of 

the metallic Impurities. The electronic property that controls the solar cell 

performance Is the thermal recombination rate of electrons and holes since 

photogenerated electron-hole pairs by solar illumination can produce electrical 

power only If they are not lost by recombination at the impurity recombination 

centers. The thermal recombination rates are represented by six parameters 

which characterize the thermal capture and emission processes of electrons and 

holes at the localized states produced by the residual impurity centers In 

a semiconductor. These are the thermal capture rates, c^ and c^, the thermal 

emission rates, e^ and e^, and the thermal activation energies, and E^^, 

of electrons (subscript n) and holes (subscript p) respectively. At thermal 

equilibrium or low applied electric field In a junction solar cell, the six kinetic 

parameters are not all Independent and only three parameters are needed. 

The others can be computed using the equilibrium electronic properties of the 

2 

semiconductor determined by other means; such as c e *c c n. where n. Is the 

n p n p 1 1 

known intrinsic carrier density; E^+E^p"Eg where is the energy gap; and 

e^“c^N^exp(-E,j^/kT) and ep*CpNyexp(-E^p/kT) where and are the effective 

density of state In the conduction and the valence band respectively. 

There are many residual impurities In a metallurgical grade silicon, such 
as A1 (7x10^^ atom/ca^), B (6x10^^), Ca (1.2x10^®), Cr (2x10^^), Cu (3x10^®), 

Fe (2x10^°), Mg (2x10^^), Mn (6x10^®), Mo (2x10^®), N1 (5x10^®), P (5x10^^), 
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Ti (2.5x10^^), V (1.2x10^^), Zr (2x10^^) and others. In addition* inpuritiea 

can b« introduced durit^ the crystal growth* such as oxygen from the molten 

silicon container made of fused quartz* carbon from the growth ambient, 

molybdenum and tungsten from the anparatus or container material and others. 

Impurities can also be Introduced during the purification of the metallurgical 

18 

grade silicon from specific purification processes such as Zn (5x10 to 
19 3 

5x10 Zn/cm )* Mo and other possible metallic contaminations from the 

crystal grower or the boat or the container for the molten silicon. Although 

many of these metallic impurities have very low distribution coefficients between 

the molten and solid silicon so that their concentration in the solid or the 

crystal can be many orders of magnitude smaller than that in the melt due 

to Impurity segregation during crystal growth. Such a reduction may still 

be insufficient to reduce the Impurity concentrations to a sufficiently 

low level so that high-efficiency solar cell can be fabricated. Estimates 

have been made for Ti based on the preliminary recombination parameter data 

which Indicated that Ti concentration must be less than about 3x10 Tl/cm'' 

in a p4’/n/nf BSF (back-surface-field) silicon solar cell in order to get 17% 

AMI efficiency. At an equilibrium segregation or distribution coefficient 

of 2xl0~^ for Ti in Si, a one-pass crystal growth or zone refining process 

19 3 

would reduce the TI concentration from 2.5 xlO Tl/cm In the feedstock Si 

to about 2.5xl0^^x2xl0"^ • 5x10^^ Ti/cm^ along 80 to 90% of the crystal leni;th 

but this is still too high to give 17% AMI efficiency. Thus, highly 

efficient and economical chemical refining processes are needed to purify 

the cetallurglcal grade silicon before it is used for single crystal growth. 

The refining processes must provide a final silicon crystal, single or poly, 

12 3 

with less than about 3x10 Tl/cro and similar amounts of other metallic liiq)urltie8 
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80 that hlgh-afflciency (17%) solar calls can be fabricated. 

To determine the electron^hole recombination parameters for each 
impurity in silicon which may degrade or control and limit the highest 
achelvable efficiency of the final solar cell requires extensive and detailed 
measurements. Because of the minute concentration of these impurities, they 
cannot be detected by conveutional analytical or mlcrochemlcal analysis and 
neutron activation techniques. In addition, these techniques would only 
give the total impurity concentration of each species and not their 
electrically active concentration nor the recombination parameters. Thus, 
alternative and electrical technlquf:s must be employed. The most sensitive 
of these are the diode transient techniques which rely on the detection of 
the minute time dependence of the concentration of the electrons and holes 
trapped at the impurity centers in the space charge layer of a semiconductor 
Junction. These include many possible variations first proposed by Sah 
and demonstrated by him and his graduate students [2] and further refined, 
extended and conq>uterized by them and other workers. A comprehensive .review 
has been given by Miller, Lang and Klmerling [3] who had themselves contributed 
to many of the refinements and automated data acquisition techniques. 

In this report, measurements of the thermal capture and emission rates 
and thermal activation energy of electrons and holes at the Ti centers in 
Si are described using two Constant Voltage Capacitance Transient (CVCT) 
methods. The results are carefully analyzed to show that they are consistent 
with the double donor model of substitutional Tl center, expected from the 
electronic shell structure of Tl atom in the tetrahedral covalent !>ond model 
of crystalline Si. 

A historical review of the recombination parameter measurement data of 
Ti center Is first given which is followed by a summary of what are believed 
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to be the most reliable or accurate results to date. A theoretical Inter- 
pretatlon of these latest results Is then given to show that they are 
consistent with the covalent bond model of silicon and the electronic shell 
structure of Tl. Here, we (.resent a new method of analyzing the theoretical charge 
states of an Impurity or defect center In a semiconductor based on the 
equilibrium statistics of the distribution of the trapped electrons and holes 
and their spin and configuration degeneracy and entropy. 



2. HISTORICAL REVIEW OF RECOMBINATION PARAMETER MEASUREMENTS 
(A) The first reported measurements of the energy levels or thermal activation 
energies of Tl related centers In silicon were given by Sah on October 12, 1976 
during the LSSA PIM (Low-Cost Silicon Solar Array Project Integration Meeting) 
sponsered by JPL (Jet Propulsion Laboratory) at Pasadena. At that time, the 
Impurity Effects on Silicon Solar Cell programs had just started with two major | ^ 

efforts undertaken by Monsanto and Westlnghouse Research Laboratories. Single 
crystal silicon Intentionally doped with metallic impurities during crystal 


growth were produced by Dow Coming for Westlnghouse and within Monsanto for 

this study. Some crystal slices were made available by JPL for our measurements. 

Two Schottky barrier runs were made quickly using a p-type silicon slice doped 

with Tl during growth from Westlnghouse (Crystal No. W-008-T1-001) . The slice 

had an elongated quarter-pie shape of approximately 1/8" by 1/2" and 14 mil 

thickness. Four-point-probe resistivity measurements gave a V/I«16.1 ohms 

for a probe spacing of 500 pm. The computed resistivity Is then p=A(V/l) 

«0.1575xl6.1>’4.06 ohm-cm which is consistent with Westlnghouse-provlded data 

15 -3 

of 4.2 to 3.1 ohm-cm. The carrier or hole concentration Is then 3.6x10 cm 

at room temperature (300^K). 30 mil diameter Mg dots were evaporated onto the 

chemically pollsed side of the surface and a blanket aluminum was evaporated 

2 

to the back side. 100x100 mil dices with four diodes per dice were cut from 


m 




T 


the wafer and mounted onto a 8-pln gold-plated TO-5 header using a low 
temperature epoxy (Epor.ek 410). A 30-mll diameter Mg dot on the top surface 
was then covered with a drop of the same epoxy and a gold preform was used to 
cover the epoxy dot. This assembly was then cured at 80*C for two hours in 
a dry nitrogen ambient. After curing, a one-mil gold wire was then bonded 
to the gold preform ultrasonically and to the post of one of the 8 pins. 

VSCTS (Voltage Stimulated Capacitance Spectra) was then taken from 77 
to 300"K which revealed three peaks at approximately 90K, 160K and 170K. The 
two lower temperature peaks were negative, revealing that they wore 
minority carrier or trapped electron emission peaks while the higher 
temperature (170K) peak was positive and due to thermal emission of holes trapped 
at a level below the midgap. The minority carrier signals and peaks were 
obtained by forward biasing the diodes which apparently gave some electron 
injection from the Mg metal into the p-Sl in the Schottky Barrier diode. 

A typical capacitance transient decay curve is shown in Fig. 1 which consisted 
of fast positive transient due to thermal emission of trapped majority carriers 
(holes in p-Sl) at the lower level and a longer negative transient due to thermal 
emission of trap/ed minority carriers (electrons) at the upper level. 

These capacitance transients were then measured at many temperatures 
around the three VSCTS peak temperatures and analyzed by a two-exponential 
nonlinear least-square-fit routine to give both the preexponential factor and 
the decay rate or thermal emission rate constants, e^ for hole emission and 
e^ for electron emission. These rate constants were then plotted as a function 
of 1000/T and least-square-fitted to the Arrhenius equations 

e^ = A (T/300)™exp[-(E„-E_)/kT] (electron emission) (1) 

n m C T 

or 

Op = Ajj^(T/300)™exp[-(E^-Ey)/kT] (hole emission) (2) 
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The experimental results were represented by the following equations 
for the thermal emission rates and activation energies* 

- 182 mV Level 

e - 3.6xlO^^(T/300)^exp(-182iaV/kT) 
n 

E^, - 470 mV Level 

e^ » 5.6xlO^^(T/300)^exp(-470mV/kT) 
ly + 336 mV Level 

ep » l.lxlO^^(T/300)^exp(-336mV/kT) 

At this eai'ller date, it was fairly certain that the two shallow levels, 

E^-182 and Gy4'336 were associated with T1 and the data were fairly accurate. 

The data for the deeper level, £^'470, were not as accurate and there were some 
doubts concerning Its origin, especially If it Is associated with surface traps 
since the Schottky barrier diode was madj on a p-type silicon surface. 

(B) The next sets of data were reported by Monsanto In their final report of 
a JPL contract, completed In Septend>er 1976 [4] but not distributed until 

late 1978. There was no capacitance transient measurements reported on their 
impurity doped silicon crystals. Impurity concentrations were determined by 
mlcroanalysls technique ? and lifetimes were determined by photoconductlve 
decay and the open-circuit voltage decay measurements. These provided some 
indications of the minority carrier capture rates at the Ti levels In the 
p-type and n-type silicon. However, the lifetimes and the capture 
rates were not and could not be correlated with the specific Tl deep levels 
since no energy level measurements were available to the Monsanto group 
before the end of the contract. 

(C) The next sets of data were reported by Northrop In their final report 
covering a contract from Decanber 1976 to October 1977 but It was not 
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distributed until Deconber 1978 [5]. In this work, the diffusion length of the 

minority carriers was measured on the Impurity doped silicon wafers from 

Monsanto and Westlnghouse/Dow-Coming, using a scanning electron microscope 

which gave very accurate, and we believe, rather reliable data. The minority 

carrier lifetimes were then computed from where L Is the diffusion 

length measured and D Is the dlffuslvlty computed from D“(kT/q)u using the 

published mobility of electrons and holes In silicon, p and p . The 

n p 

diffusion-length lifetime data were also correlated with the photoconductlve 
decay and open-circuit voltage deca. lifetimes obtained by Northrop as well as 
reported by Monsanto and tfestlnghouse on their Impurlty-doped silicon crystals. 
This correlation curve was the basis for calculating the minority carrier 
capture rates at the T1 levels to be presented In subsection E by us. Again, 
qo correlation of the lifetime data were made with the specific T1 levels since 
reliable and accurate Ti level Information were not available at that time. 

(D) The next set of data on the thermal activation energies as well as the 
thermal emission and capture rates of electrons and holes were reported by 
Westlnghouse in their 11-th Quarterly Contract Report in July, 1978 [11] . 
Measurements were made for many impurities using the VSCTS (Voltage Stimulated 
Capacitance Transient Spectroscopy) . It Is also known as DLTS - Deep Level 
Transient Spectroscopy or JTS - Junction Transient Spectroscopy but DLTS and JTS 
are all Inclusive names for all the possible variations first reported by Sah 
and his graduate students [2] and reviewed by Miller, Lang and Klmerllng [3], 
while VSCTS Is a specific member of the DLTS or JTS family. There were 
considerable variations among these first atempts as revealed by the table 
given on the next page. Part of these were probably due to the uncertainty in 
the identification of the VSCTS peaks and to noise and calibration accuracy of 
the early Instrumentation. These were mostly resolved In their later data 
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reported in their 16-th Quarterly Contract Report to be described in subsection 
(H) and in their papers to be described in subsections (F), (G) and (I). 

TABLE 1 RECOMBINATION PROPERTIES AT T1 CENTERS IN Si 
(Westlnghouse Data, July 1978, p.llO (6J) 


Crystal 

No. 

Diode 

Type 

Energy Level 
(mV) 

Density 

(cm“^) 

finlsslon 

Cross-Section 

(cm^) 

Emission 

Rate 

(cm^/sec) 

W-008 

m/p 

Ec - 340 

7.0E14 

4.0E-16 

4.0E-9 



Ey + 290 

7.0E13 

4.7E-14 

4.7E-7 

W-008 

n/p 

Ec - 260 

1.5E13 





Ey + 290 

2.6E13 

? 


W-033 

?/p 

Ec - 270 

4.0E12 

4.5E-17 

4.5E-10 

W-065 

?/n 

Ey + 450 

3.6E13 

3.6E-16 

3.6E-09 



Ey + 250 

3.0E13 

1.2E-15 

1.2E-08 


(E) A comprehensive analysis was then made by Sah of all the published data 
to date and reported in March, 1979 [7]. At that time, the most accurate thermal 
emission and activation energy measurements were those reported by Sah and 
described in subsection (A). In addition, a detailed analysis was also made of 
all the minority carrier lifetime data reported by Westlnghouse, Monsanto and 
Northrop. The lifetimes, corrected using the correlation curve between Che 
diffusion length lifetime and the photoconductive and open-circuit voltage decay 
lifetimes were plotted as a function of the T1 concentration, either measured 
by the microanalysis techniques or estimated from an assumed liquid-solid 
segregation coefficient. The t vs N^^ data of twelve Ti-doped n-Sl and p-Si 
from Westlnghouse and four from Monsanto were plotted on the logr vs logN^^ scale 
in Fig. 2. Since r-CcN^^)”^ or logT*-logC-logN^j^ where C is the thermal capture 


- 93 - 


LIFETIME (sec) 



TITANIUM CONCENTRATION (Ti/an3) 


Fig. 2 The minority carrier lifetime computed from diffusion length 
measurements as a function of the Ti concentration from 
mlcroanalysls of many p-type and n-type Si crystals doped with 
Ti during growth. Solid dots are p-Si and electron lifetime and 
circles are n-Si and hole lifetime. Numbers with character M or 
C in front are crystals grown by Monsanto while those without are 
grown by Dow-Corning/Westinghouse. 
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rate of the minority carriers, three heat-gueae llnea were drawn In Fig. 2. Three 


lines were chosen because the points seem to group along the three lines shown 


In Fig. 2 and because the preliminary measurements gave three thermal activation 


energies as described in subsection (A). Later analysis indicated that the 


deep level at E^-470 mV was due to a thermally Induced center not related 


directly to the double donor Ti center. From these three lines, three 


minority carrier capture rates could be calculcated at the lifetime measurements 


temperature (300”K). The results [7] are summarized below. 


E^-182 Acceptor Level 


c * 1 . 3x10 
n 


(cm^/s) 


e - 3.6xl0^^(T/300)^exp(-182/kT) - 3x10® ( 1/s) 

n 


Cp • 5. xio ’ (Fig. 2) 


(cm®/s) 


e - 2.2x10"^ (297 *K) 

P 


( l/s) 


E--470 Donor Level 

c 


c *2* xlO 
n 


e » 5.6xlO^®(T/300)^exp(-470/kT) 
n 


Cp - 2.5x10 ° (Fig. 2) 


e » 9.2 
P 


(297*K) 


Ey+336 


Donor Level 


c » 3. 7x10' ' (Fig. 2) 
n 


e “1.6 
n 


(297“K) 


c » 8.2x10 ’ 

P 

e - l.lxl0^^(T/300)^exp(-336/kT) » 2x10^ 
P 


The numerical values shown above were computed at 297*K. The two upper levels. 


E -182 tnd E„-470 mV detected from minority carrier emission transients were 

c c 
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tentatively attributed also to the T1 center with the shallower level » Ec-iS2. 

designated as an acceptor while the deep level, Ep*>470, as a donor based on 

2 2 

the electronic shell configuration of Tl, (Ar)3d 4s , whose two 4s electrons 
are weakly bounded and can give rise to the two donor states at E^-470 and Ey-f336. 
It was further speculated that the shallow E^-182 level may be the oxygen- 
vacancy level rather than Tl related. Later results Indicate that the 
Identification of the E^+336 level as the second Tl donor level Is correct but 
the first Tl donor level lies at E^-228.65 mV instead of E^-470 mV so that 

the level at E.-182 mV was probably the second Tl donor level. The difference 

c 

In energy (182 vs 228.65) was partially due to the Inaccuracy of the original 
data analysis since It was a minority carrier transient and might have some field 
dependences. The data were used to compute the characteristics of Tl-doped 
solar cells made from the 14 Westlnghouse crystals (7]. The measured and 

computed J._, V*-, EFF were In good agreement but the agreement with the 

SC oc 

measured fill factor was poor. It was speculated that this poor agreement 

was probably due to a la.'ge contact resistance which was not taken Into account 

In the theoretical transmission line equivalent circuit model. It was also 

noted that the agreement was good at low Ti concentrations but poor when the 

Ti concentration reached 4x10 Tl/cm in crystal W-008. It was suggested [7] 

that Impurity segregation and cluster may have reduced the experimental 

J anH EFF. It was also noted this was the first Ti-doped Si crystal grown 
SC 

and its quality may have been affected by some unknown and uncontrolled growth 
parameters. In any case, the computer model was demonstrated to give good 
correlation with actual experimental solar cell data as a function of the Ti 
concentration, even though the recombination model was still uncertain, 
because the recombination parameters were accurately measured and used as the 
Input parmneter for the computer calculations. It also demonstrated that to get 
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a 17Z AMI efficiency cell, the T1 concentration auat be leea than 2.5x10 -cm 
and a p+/n/n^ cell structure must be used since the n^/p/p^ cell structure was 
less efficient because electron lifetime In Tl-doped p-Sl Is lower than hole 
lifetime In Tl-doped n-Si. 

(F) The recombination properties of the Ti centers In SI were extensive 

and carefully Investigated by many experiments carried out by Chen, et. al. 

a« Camegle-Mellon and reported In November, 1978 which was published In 

September 1979 [8]. Both Westinghouse Tl-doped SI and Ti diffused SI were 

used to fabricate both diffused and Schottky barrier diodes on both n»Si 

and p-Sl. Three electron and three hole trapping levels were observed with a 

40 MHz VSCT Spectrometer using double boxcar correlation technique. Two of 

these six levels were shown to be directly proportional to the TI concentration 

and were attributed to a single TI center. They are E^-264 mV and E^4-290-35 mV. 

They were both designated as donor levels since a 360 ohm-cm n~Sl was 

diffused by Ti and did not show acceptor or p-type conversion. The other four 

levels were attributed to thermal generation at high diffusion temperature or 

unintentionally incorporated impurities during processing. 

Detailed majority carrier capture rates were also measured which 

Indicated that hole capture into the second Ti donor level, Ey-f-290 mV was 

thermally activated with an activation of 35 mV. This gave an effective 

thermal activation energy of holes of Ey+29i>-35*Ey+255 mV. 

Hall effect measurements were also carried out for the n-Sl sample 

giving a electron activation energy of E^-220 mV Instead of the E^-264 mV 

determined by VSCTS method. The capture rates and cross sections at 300*K are 
—IS 2 — fl 3 

0 “3. 5x10 ^ cm and c *3x10 cm /s for the first donor level at E_-264 mV 

nc n C 

—17 2 —10 3 

and a^^*l«7xl0 cm aad c^*2xl0 cm /s for the second donor level at mV* 
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UDfortuTiateJy, detailed leeet-equares-flt reeulte of the emlaslon rate data 
were not ,qlven and only a graphic cosq>arlson could be aade with earlier and 
subeequent and more accurate data. 

(G) More detailed capacitance transient measurements were made by Westlnghouee 
and these were reported In three publications, a paper In the Solid State 
Electronics Journal [9], their 17-th Quarterly Contract Report [10], and a 

IEEE Transaction on Electron Devices Journal article [11]. The two levels 
identified by Chen as those from T1 are now reported also by Westingbouse 

who also gave values for the majority carrier capture rates In two of these 

•18 2 •!! 1 

references. These were Op^-5.3xl0 cm or Cp«5.3xl0 cm /s at the £^+300 mV 

I/O TO 

level [10] and 0 -1. 04x10” or c ■1.04x10” ^cm'*/e at the E -260 mV level and 
nc n C 

•16 2 —9 3 

Op^«8.3xl0 cm or Cp»8*3xl0 cm /s at the Ey+300 mV level. These values were 

•18 ? —16 2 

self-inconsistent (for example, a,^*5.3xl0 cm' and 8.3x10 cm or a factor 

pc 

156.6 dlfferoat) and they were different from those given by Chen which was 
described In sub-section (E). The values given by Chen are considered more 
reliable. 

(H) The T1 levels were also measured by a JPL in-house group [12] using 
solar cells fabricated by Westlnghouse on the Dow-Cornlng Tl-doped p-type Si. 
Measurements were made on either mesa or photo lithographically etched nf/p 
diodes. The VSCTS data showed three majority carrier peaks and one minority 
carrier peak, at Ey+0.18, Ey+0.27, Ey+0.37 and E^,-0.24 eV. Only the Ey+0.27 
level appears In the original wafer while the other appears after high 
teaq>erature diode or cell diffusion. The other levels were attributed to 
Titanium-Oxygen complexes since titanium Is highly electronegative. This model 
Is similar to the Alu m i n u m -Oxygen complex which was observed In A1 -doped 
silicon [13]. 


(I) More extenelve cepeeltaace trenelent epeccroecopy wee also obtelned (14) by 

Cheng and Leung In a J?L In^^ouse work uaing the Weetinghouee Tl-doped p-Sl 

aateriala and eolar cella aa well as silicon ribbons, by the web dendritic 

process (Wsstinghouse) and edge-defined-growth process. The results are 

o^"5xl0”^^ cm”^ and Op«5.5xl0’’^®cm”^ for the E^+280 mV T1 level at 200K. 

-15 2 

They also reported an estimate for the E.-240 mV of o^ > 4x10 cm . 

(J) We have carried out extensive and highly accurate measurements of the 
thermal emission rates and activation energies of trapped majority carriers 
at the T1 levels in n-type and p-type Si prepared by Westinghouse/Dow-Corning 
as well as Ti-doped Si diffused and Schottky Barrier diodes provided by 
Westinghouse. These results are sianmarized below. 

+ 305 mV 

e - 5.45xl0^^(T/300)^exp(-3O5.34/kTj (1/sec) (1) 

P 

\dtere the reduced chi square is 1.036 and the rms deviation of the two constants 
are 5.45ll«13 and 305.34±2.31« The large uncertainty in the pre-exponential 
factor was due to the msall nuifter of temperatures (nine). The results are 
accurate since the rms deviation of the emission rate is only 1.99Z. 

E,, - 231 mV 

e - 1.42xl0^°(T/300)^exp(-231.34/kT] (1/sec) (2) 

n 

%diere the reduced chi square is 0.73 and the rms deviation of the two constants 
are 1.4210.080 and 231.3410.54. The results are quite accurate and has a 
rms deviation of only 1.16Z. The capacitance decay curves were taken at 
thirteen temperatures in the range of 101.60 to 121. 81K. 

The thermal capture rate measurements have not been obtained for all 
of the electron and hole transitions. The majority carrier capture rates can be 
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obtAlMd by fMrclal»tr«p-filllQg «q>erinants over a wide rente of teoq>erecuree 
below Che aMucimiB tanperacure where Che chemal emlaaion beconee coo feac Co 
be aeaaured by Che capacicaoce aecer. Above chla cemperacure, current 
cranalenc imiac be nonlcored under parclal-crap-flllir^. The minority carrier 
capture ratea were meaaurml either by Injecting mlnorlcy cnrrlcrs Into the epace 
charge region ao chat they can be captured to give a Junction transient or by 
Junction current switching transient. Thus far, measurements are siade on the 
hole capture rate at the second donor level* Ey‘f305 mV* using the technique of 
partlal*'trap*‘fllling by majority carriers. The data of our experiments as well 
as those of others are given In Table 2 and plotted In Fig. 3. 

It Is evident that there are wide disagreements among the hole capture 
rate data glvra by Weatlnghouse's tvo reports and Carnegle-Mellon as noted 
earlier. Our data at much lower temperatures than those of Carnegle-Mellon 
appear to lie on a reasonably smoorJi curve In Fig. 3. Tliere are some 
spreads In our data taken on different diodes which are illustrated in Fig. 3 
where the data symbols are indicated in Table 2. are for the diffused nf/p 
diode TiP152Dl;^are for the Schottky barrier diode T1P137S1; and ^ are 
for the Schottky barrier diode T1P152S1. The following results are evident. 

(1) c is the lowest in the Schottky barrier diode, 71P137S1 ( ^ ) which has 

P ^ 

the highest boron concentration among these three diodes (4.69x10 boron/ cm ). 

(2) Cp is higher in the SB diode T1P152S1 (^) tl»an in the diffused n-f/p 
diode T1P152D1 (@ ). (3) All three sets of data of c^ from our measurements 
are higher than the values extrapolated from the Carnegie^llon-Westinghouse 
high- tempera Cure data ( • ). 

One important point to note Is that the thermal activation energies 
and emission rates of our three diodes were all nearly equal to each other* 
assuring that the capture measurements were at the same Ti center. 
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Table 2 Hole Capture Rate at T1 Center in p-Si 


Diode Ho. 

Temperature 

1000/T 

3 

c (cm /sec) 

Source of Data - Conments 


(K) 


-pj 

This work. ® In Fig. 3. 

T1P152D1 

77.01 

12.98 

1.620E-11 

(n+/p) 

90.00 

99.66 

11.11 

10.03 

2.477E-11 

3.018E-11 

N^"5. 96x10^^ cm~^ 


109.99 

9.092 

3.608E-11 



128.92 

7.757 

4.733E-11 



136.95 

7.302 

5.373E-11 



144.95 

6.900 

5. 5052-11 

This work. A in Fig. 3. 

T1P137S1 

109.82 

9.106 

2.532E-11 

(m/p) 

129.05 

136.75 

7.749 

7.313 

3.617E-11 

4.434E-11 

N^“4. 69x10^^ cm ^ 


144.90 

6.901 

5.454E-11 

Large series resistance 
prevented low temperature 
measurements . 

T1P152S1 

77.01 

12.98 

1.931E-11 

This work, in Fig. 3. 


109.87 

9.102 

4.593E-11 



144.98 

6.898 

6.010E-11 


Diffused & 

179.89 

5.559 

0.761E-10 

Chen-Mllnes-Rohatgl [ 8 ] . 

Schottky 

196.14 

5.098 

0.957E-10 

# In Fig. 3. Read off the 


201.90 

4.953 

1.03 E-10 

figure from [8J as accurately 


213.80 

4.677 

1.17 E-10 

as possible. 


221.43 

4.516 

1.29 E-10 



241.66 

4.138 

1.54 E-10 



250.25 

3.996 

1.68 E-10 



2 


8.3 E-09 

Westlnghouse [11]. 


2 


5.29 E-11 

Westlnghouse [10]. 


200 

5.000 

5.05 E-11 

JPL Cheng-Leung [14]. 


CAPTURE RATE (cmVs) 


ORIGINAL PAGE IS 
OF POOR QUALITY 



RECIPROCAL TEMPERATURE 
1000/T (k’’) 


Fig. 3 Temperature dependence of the hole capture rate at the 
second Ti donor level, E^+305 mV, in p-type silicon. 
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The capture rate plotted In Fig. 3 la the quantity measured directly 

from the experiment. The data presented by Chen-Mllnes-RohaCgl In the 

Camegle~Mellon paper [8] was given as capture cross section. A 

conversion back to the capture rate required the use of an thermal velocity 

expression for holes since It is temperature dependent. The rms velocity, 

0p«/3kT/mp, is used where nip Is the density-of-state effective mass of the 

valence band . nip Is slightly tonperature dependent and Is ctmiputed using 

nip/m « 0.64742[1+(T/193.472)3 ]x[1+(T/142.739)3]/{[1+(t/202. 056)3] [1+(t/166. 785)3]}. 

The capture rate is computed from the measured filling time constant using 

c^*(Pt^) ^ where Is the filling time constant and F is the hole concentration. 

Deionization correction is also made to compute F from the boron concentration, 

_2 

^AA measured from C vs V at 300K. The formula for computing F is 

F - where K^«(N^/g^)exp[-(E^-E^)/kT] , and 

N^“2. 51x10^® (T/300)^^^(nip/m)^^^. E^-E^=0.0444 eV and g^-2 are assumed. The 
deionization correction is significant but not very large. 

If a line is drawn through these data, a rough or poor fit can be made 
to give a Arrhenius equation for the hole capture cross section which is 

- 6.2x10'^^ exp(-34.85/kT] cm^ (3) 

This would change the thermal activation energy for hole emission by 34.85 mV. 

The electron capture rate has not been measured because it is very 
large and the filling time required is below 10 ns making it impossible to 
perform the partial filling experiment due to pulse shape distortion by the 
loading of the Tl-doped diode. This filling time can be increased substantially 
if the resistivity of the Tl-doped n-type Si crystal is Increased. At 30 ohm-cm, 
the hole concentration will be sufficiently low to allow for accurate measuronents 
but high resistivity Tl-doped n-type silicon crystal is not available.^ 

# It became available after the completion of the draft of this technical report. 
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TBBORSTICAL ANALYSIS 


Although Che complete set of four capture rates at the two T1 levels 
have not been obtained, the data on hand for the majority carrier capture rates 
can be used to delineate the electronic models and the charge states of the Ti 
levels In p-type and n-type silicon. This Is possible by noting that at thermal 
equlllbrlisB, the emission and capture rates are related by the mass action law. 
For holes crapped at the lower Tl level, this Is 

e « c N^expl-(E^-Ey)/kT] - a 0 IJ^exp(-(Ej-F^)/kTl (4) 


where we have used the definitions for the cross sections given by 


c » 00. (7) 

n n n 

Here, a =/3kT/m and a ■•/SkT/m are the rms velocity of the holes and electrons 
p p n n 

respectively whose effective masses are m^ and m^. In our previous data conversion, 
we have employed the denslty-pf-state effective masses, sip and m^. 

In the experiments of majority carrier capture and emission, such as 
holes at the lower Tl level, Ey+SOS mV, both e^ and c^ can be measured Independently 
but they are measured under drastically different conditions. The emission 
transient and rate Is due to the thermal emission of the trapped carriers from 
the Tl centers In the junction space-charge region where there Is a high electric 
field and hence It Is a highly non-equlllbrlum process. The capture transient 
and rate Is due to the thermal capture of carrier by the Tl centers In the 
quasl-neutral region and hence It Is an equilibrium process. Because emission 
Is not at thermal equilibrium, we define a separate emission cross section to 
distinguish It from the capture cross section. Thus, the four equations for the 
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four experimental meaauremente are 


and 


•p ■ «p.»p'V“i't-(v®v)/“i 


c® • a* 0 

P pc p 


•n ■ W'”' 


c* - o* 9 . 

n nc n 


( 8 ) 

(9) 

( 10 ) 
( 11 ) 


which give four experimentally determinable cross sections, and for 

holes and and o®^ for electrons where the superscript ® denotes equilibrium. 

In addition to the equllibrlum-nonequllibrlum difference, the two 

cross sections for holes (or electrons) are also related through a factor vdiich 

depends on the nature of the charge state or impurity potential which determines 

the spectrum of excited bound states, their spin and configuration degeneracy, 

as well as a conflguratlcn entropy due to lattice relaxation around the impurity. 

We shall develop or summarize these theoretical results which relates <7^^ to 

0 and a to a . Thus, the experimental data can be used to determine which 
pc ne nc 

theoretical model is consistent with the measurements which would then allow a 
determination of the charge state or the Impurity potential. 

We shall first obtain the complete expression for the emission rate 
at thermal equilibrium in terms of the thermal equilibrium capture rate or cross 
section for a donor impurity level. Suppose that the donor impurity has five 
valence electrons with a ground state energy level at and a degeneracy 

of gp. Then, from the equilibrium mass action law or from balancing the 
electron emission and capture transition rates, ®n**D“c*NPjj=c^N(Npjj-Njj) and 
using the equilibrium distribution function for the trapped electrons [15] 

Up - Njjp/{l+g"^expI-(Ep-E^)/kT]} (12) 

and the electron concentration given by 

N - N^exp[-(E^,-Ep)/kT] (13) 

- 105 - 


we have 


Comparing the non<-equlllbrltBn experimental relationship, (10), with the 

theoretical equilibrium relationship, (14), we can write 

- g /r (16) 

nc ne D n n D n 

where r =e /e* > 1 since the emission rate In a high field will be larger due to 
n n n — 

binding potential barrier lowering and tunneling. For this donor impurity, 
the hole emission and capture cross sections are related by 




(17) 


idiere r =e /e* > 1. 

p p p- 


A similar derivation can be made for an acceptor Impurity level whose 
ground state has a g^-fold degeneracy. The results for the acceptor impurity 
are different than those for the donor loqiurlty given by (16) and (17) since 
holes are trapped at the acceptor level while electrons are trapped at the donor 
level. The cross section ratios at an acceptor level are given by 


and 


o® /a ■ g./r 
pc pe A p 


o® /a “ (g r )~^. 
nc ne A n 


(18) 


(19) 


If we take Into account of the excited states and their degeneracy, then 
the degeneracy factor, g^, must be replaced by a sum-over-state [15,16] given by 

*D ■ «1>1 * + — • «« 

where g|^^ and (1«1,2,3, . . . ) are the degeneracy and energy level of the 1-th 
excited state. The ground state Is !■!. 
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The different depmdences of the capture-remission cross section ratio 


of a donor and an acceptor Impurity* given by (16) and (17) for a donor and 

(18) and (19) for a acceptor* makes It possible to distinguish them from 

experimental measurements of the cross sections from the capacitance transient 

experiments. To demonstrate this possibility* we tcdice the simplest situation 

In which there Is no electric field dependence of the emission rate so that 

r *r ■!. Then* we have 
n p 


Donor Level (Electron Trap) 

i «D - 2 


a* /cr 
nc ne 




Acceptor Level (Hole Trap) 

/o - )”^ < g"^ < 1/2 

nc ne _A n — A - 


0® /a 
pc pe 




^8. i 2 


( 21 ) 

( 22 ) 

(23) 

(24) 


These sets of equations show that the experimentally determined ratio of the 

capture to emission cross section will make It possible to determine If the level 

Is a donor level or an acceptor level. 

Let us apply 'the above result to the measurements made on the T1 levels. 

From the hole emission rate at E^i*305 mV given by (1), the hole emission 

-15 2 

cross section at 300®K is a *2.67x10 cm . The hole capture cross section 

pe 

—17 2 

calculated from (3) Is o® *1.6x10 cm . The ratio Is o® /o *1/167*0.006. 

pc pc pe 

Thus, this lower T1 level must be a donor level. The small ratio Indicates that 
r *e /e® » 1 since g_ in (22) cannot be substantially greater than about 2. 

p p p “T) 

Part of this large ratio is due to electric field dependence but configuration 
entropy due to strong electron-phonon coupling and lattice distortion Is likely 
to contribute. If It Is due entirely to the latter* the entropy Is klog_166BSk 
which Is not unreasonable. 
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A similar calculation may bo mada for the electron emission rate at 

Ep~231 mV glvm by (2). The emission cross section at 300*K Is a *4. 46x10*^ ^cm^. 
u ne 

The electron capture rate and cross section have not been measured In the 
Tl**doped n-Si but It has been measured In n-Sl where T1 was Introduced by 
diffusion. This was reported by Chen<^lnes-Sohatgl (8] to be a constant 
from 150 to 210*K with a value of o*^«3.5xl0”^^ cm^. Using this value* the 
ratio is then with (21) to (24) shows that this 

T1 level must also be a donor level. 

Thus, from the emission and capture cross section measurements of 
majority carriers at the two T1 level* we can conclude without doubt that the 
two T1 levels are donors. The most likely possibility is that T1 is a double 
donor with cliarge states 0* ■t'l and +2 and the upper level* E.-231 mV the 
first donor level and the lower level* E^*f305 mV the second donor level. 

The large ratios of the capture to emission cross sections suggest 

that there is considerable lattice relaxation from strong electron-phonon 

interaction since electron spin and electron configuration degeneracy cannot 

account for such large ratios neither can electric field enhancement of the 

hole emission rate at the lower level. Electric field enhancement of the 

electron emission rate at the upper level would give the wrong change since 

it would increase r or e /e^ and make the a* /a ratio even smaller, 
n n n nc ne 

The conclusion here that the two T1 levels are donor-1 Ike are 
consistent with the conductivity-type conversion experiment of Chen-Mllnes-Rohatgi 
in which they find no conversion of a 320 ohm-cm n-Sl diffused with Tl* 
implying that all Ti levels must be donor like. The recombination coefficient 
measurements given here not only determines whether Tl Is donor or acceptor but 
also provide additional indication that strong electron-phonon interaction 
occurs around the Tl impur:ity site to give a large configuration entropy. 
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The conclusion about the double donor nature of the T1 impurity in Si 

gives further support to the valence bond model. The electronic shell structure 
2 2 

of Ti lb [(Ar)3d 4s ] The radii of the tvo 4s electrons are large, 1.477A, 
compared with one-half of the nearest neighbor silicon spacing (2.3S/2«1.17A) . 
Thus, these two 4s electrons of the Tl atom are loosely bound in the Si 
crystal lattice and they are easily released to produce two bound states, the 
E^-231 and the E^-t‘305 mV levels. 

The electron capture cross sections and rates at these two donor levels 
are expected to be large due to attractive Coulomb potential. On the 
other hand, hole capture cross sections are expected to be small due to the 
repulsive Coulomb potential. Thus, the electron lifetime in p-Sl should be 
considerably shorter than hole lifetime in n-Sl, making the p4*/n/nrf silicon 
solar cell containing 7.'1 likely to have higher efficiency than the n+/p/p+ 

This was verified by computer design simulation [7]. 

A similar conclusion on solar cell efficiency can also be reached for 
the vanadium recombination center which is a double donor in silicon and other 
double donor recombination centers whose two levels are located in the two 
different halves of the silicon energy gap. 
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